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ABSTRACT 

The explosion of ultra-stripped stars in close binaries can lead to ejecta masses <0.1 Mq 
and may explain some of the recent discoveries of weak and fast optical transients. In Tau¬ 
ris et al. (2013), it was demonstrated that helium star companions to neutron stars (NSs) may 
experience mass transfer and evolve into naked ~ 1.5 Mq metal cores, barely above the Chan¬ 
drasekhar mass limit. Here we elaborate on this work and present a systematic investigation of 
the progenitor evolution leading to ultra-stripped supernovae (SNe). In particular, we examine 
the binary parameter space leading to electron-capture (EC SNe) and iron core-collapse SNe 
(Fe CCSNe), respectively, and determine the amount of helium ejected with applications to 
their observational classification as Type lb or Type Ic. We mainly evolve systems where the 
SN progenitors are helium star donors of initial mass Mn e = 2.5 — 3.5 Mq in tight binaries 
with orbital periods of P or b = 0.06 — 2.0 days, and hosting an accreting NS, but we also 
discuss the evolution of wider systems and of both more massive and lighter - as well as sin¬ 
gle - helium stars. In some cases we are able to follow the evolution until the onset of silicon 
burning, just a few days prior to the SN explosion. We find that ultra-stripped SNe are possible 
for both EC SNe and Fe CCSNe. EC SNe only occur for Mne = 2.60 — 2.95 Mq depend¬ 
ing on P 0 rb- The general outcome, however, is an Fe CCSN above this mass interval and an 
ONeMg or CO white dwarf for smaller masses. For the exploding stars, the amount of helium 
ejected is correlated with P or b - the tightest systems even having donors being stripped down 
to envelopes of less than 0.01 Mq. We estimate the rise time of ultra-stripped SNe to be in 
the range 12 hr — 8 days, and light curve decay times between 1 and 50 days. A number of 
fitting formulae for our models are provided with applications to population synthesis. Ultra- 
stripped SNe may produce NSs in the mass range 1.10 — 1.80 Mq and are highly relevant for 
LIGO/VIRGO since most (possibly all) merging double NS systems have evolved through this 
phase. Finally, we discuss the low-velocity kicks which might be imparted on these resulting 
NSs at birth. 

Key words: supernovae: general — binaries: close — X-rays: binaries — stars: mass-loss — 
stars: neutron — white dwarfs 


1 INTRODUCTION 

The nature of a supernova (SN) in a close binary determines if 
the system will remain bound or dissolve into two runaway com¬ 
ponents after the explosion. The two main decisive factors are: 
i) the kick velocity added to the newborn neutron star (NS) rel¬ 
ative to the pre-SN orbital velocity of its progenitor star, and ii) 
the amount of mass ejected during the SN. For bound systems the 
post-SN orbital characteristics and systemic velocities depend on 
these quantities which are therefore related to whether the explo¬ 
sion is an iron core-collapse SN (Fe CCSN) or an electron-capture 
SN (EC SN). It has been demonstrated that Fe CCS Ne may p ro- 
duce kicks up to, and even above, 1000 kms -1 ( IScheck et al .1 
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l200d : ljankall2012l : Wong wathanarat et~aill20 1 3l) whereas EC SNe 
are expected to be almost symmetric (i.e. fast explosions) and pro- 
duce very small kicks of less than, perhaps, a few 10 kms -1 (e.g . 

IPodsiadlowski et all2004l : lKitaura et all2006l : lDessart et al.l2006f) . 

It is well known that the outcome of stellar evolution in close 
binaries differs significantly from that of single stars. The main ef¬ 
fects, besides from mass loss/gain, are changes in the stellar rota- 
tion rate, the nuclear burning scheme and the wind mass-loss rate 
dde Minkll2O10l: lLangerll2012) . As a result, bin ary interacti ons af- 
fect the final core mass p rior to the collapse dBrown et alJl200ll : 
IPodsiadlowski et alj2004 ), and therefore possibly the nature of the 
SN event, the type of compact remnant left behind, the amount of 
envelope mass ejected, and thus the kinematics of the post-SN bi¬ 
nary. Given that almost all massive stars are members of close bi¬ 
naries dChini et alj2012l) and 70 per cent of them interact with their 
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companion star l lSana et al.ll20 1 2j| . it is important to probe the evo¬ 
lution leading to SNe in close binaries. 

In recent years, high-cadence surveys and dedicated SN 
searches have increased the discovery rate of unusual optical 
transients with low peak luminosities (< 10 42 ergs -1 ) and/or 
rapidly decaying l ight curves. Exa mples comprise SN 2 005ek 


(Drout et al 


20131). SN 201 OX l lKasliwal et al.1 l201dh and 


SN 2005E i Perets et ai]|2010l) . These SNe show diverse spectro¬ 
scopic signatures (cf. Kleiser & Kasenll2014h and thus their expla¬ 


nation may require a diversity of models. iDrout et al.l J 2 OI 3 I ) con¬ 
cluded that SN 2005ek represents the smallest ratio of ejec ta to 
remnant mass observed for a core-collapse SN to date. lTauris et al .1 
d2013l) recently suggested that ultra-stripped SNe may be respon¬ 
sible for such events and demonstrated that synthetic light curves 
of these SNe, in which only ~ 0.1 Mg of material is ejected from 
an exploding star barely above the Chandrasekhar mass limit, are 
consistent with observations of SN 2005ek-like events. 

Furthermore, merging NS/black-hole binaries are the prime 
candidate sources of high frequency gravitational waves to be 
detected by the LIGO/VIRGO network within the next 2-3 yr 
dAasi et all2013l) . In th is connection, a search for th eir electromag¬ 
netic counterparts (e.g. iKasliwal & Nissankell20l4 ) is of great in¬ 
terest. A good understanding of the formation of such systems, the 
nature of the SN explosion which created the last formed NS and 
its expected light curve, as well as the distribution of the resulting 
NS masses, is therefore of uttermost importance. 

Most Type Ib/Ic SNe are expected to originate in bi¬ 
nary systems from the initially more massive star which has 
lost its hydrogen envelope by mass-transfer to its companion 
dPodsiadlowski et all 1 19921 : lEldridge et aid |2008|) ; th ese pre-SN 
stars typically have an envelope mass of 1 Mg or more dYoon et al .1 
l20lrf) . In close X-ray binaries, however, a second mass-transfer 
stage from a h elium star to a NS can strip the helium star further 
prior to the SN dHabetslll986l:lDewi et alj2002l:lDewi & Polsl200i 


llvanova et all2003l : lTauris et all201 3 L and references therein). The 

fast decay of the Type Ic SN 19 941 with an ab s olute magnitude of 
My « —18 was explained by iNomoto et aid dl994l) : Isauer et al.l 
d2006l) from such a model, for which they suggested a carbon- 
oxygen star progenitor with a mass of ~ 2 Mg and a correspond¬ 
ing ejecta mass of ~ 0.9 Mg. As mention ed above, an even more 
extreme case was demonstrated by iTauris et~aT] d201 3). hereafter 
Paper I, in which only ^0.1 Mq of material was ejected. A simi- 
larly small ej ecta mass was previously sugge sted by lPiran & Shavfvl 
d 20041.120050 and lPodsiadlowski et all 120051 ) for the second formed 
pulsar in PSR J0737-3039. __ 

A common envelope (CE) phase (e.g. llvanova et alJl2013h re¬ 
sulting from high-mass X-ray binaries (HMXBs) is thought to have 
formed the tight helium star-NS binaries, which are the starting 
point for the above-mentioned calculations, leading to a second SN 
explosion in these systems. For a general rev iew on th e forma¬ 
t ion a nd evolution of HMXBs, see e.g. ITauris & van den Heuvell 
d2006h . Of particular interest for the scenario discussed in this 
paper is mass transfer by Roche-lobe overflow (RLO) initiated 
by the helium-star expansion after cor e helium depletion, i.e. so - 
called Case BB RLO, first dis cussed bvISavoniie & Taken j I 1 9761 ): 
iDe Greve & De Loorel d 19771) : lDelgado& Thomasl 1 1981 1. 

Here we investigate the progenitor evolution in such close bi¬ 
naries leading to a SN event. The aim is to determine which stellar 
systems lead to EC SNe, which systems produce Fe CCSN. and 
which systems avoid a SN altogether and leave behind a massive 
white dwarf (WP) remnant. We accomplish this by a systematic 
and detailed modelling of about 70 helium star-NS binaries in close 


orbits. In Section [2] we briefly refer to the stellar evolution code 
behind our numerical computations, and in Section [3] we present 
our models with detailed examples of four different outcomes. We 
highlight a peculiar model in which extensive mass loss results in 
the complete quenching of the central nuclear burning. In Section[4] 
we discuss the expected observational properties of our exploding 
stars and in Section [5] we discuss the mass-transfer rates from the 
helium star donors. In Section[6] we discuss the importance of ultra- 
stripped SNe for the formation of double NS systems and comment 
on their masses and imparted kicks. In addition, we estimate their 
merger times which are relevant for LIGO/VIRGO. Further discus¬ 
sions of our findings are given in Section[7] including a comparison 
with previous work and a mapping of initial parameter space to fi¬ 
nal outcome. We summarize our conclusions in Section [8] and in 
the Appendix we analyse the timescale for NS in-spiral in the dilute 
envelope of a pre-SN star. 


2 BINARY STELLAR EVOLUTION CODE AND INITIAL 
PARAMETER SPACE 

For our detailed calculations of the evolution of helium star- 
NS binaries, we applied the binary stellar evolution code BEC 
(e.g. lYoon et al J 1201 Ol . and references t herein ). This code was 
originally develope d by Wellstein et all J200 ll) on the basis of 
a single-star code dLangeill 1 9981 . and references therein). It is a 
one-dimensional implicit Lagrangian code which solves the hy- 
drodynamic form of the stellar structure and evoluti on equations 
dKippenhahn & Weigertlll99(il : iKozvreva et alj f2014h . The evolu¬ 
tion of the donor star, the mass-transfer rate, and the orbital sep¬ 
aration are comp uted simultaneously throug h an implicit coupling 
scheme (see also IWellstein & L ange ill 1999 ) usin g the Roche ap¬ 
proximation in the formulation of Eggletonl ' 1983^. To com pute the 


mass-transfer rate, we use the prescription oflRitten 1 9881) . We em¬ 
ploy the radiative opacities of llglesias & Rogers! ( jl99fih . which we 
interpolated in tables as function of density, temperature, and chem¬ 
ical element mass fractions, including c arbon and oxygen. For the 
electron conduction opacity, we follow [H ubbard & Lampel i ll 9691) 
in the non-relativistic case, and ICanutd dl97ol7 in the relativistic 
case. Th e stellar models are computed using extended nuclear net¬ 
works dHeger et ahlboooi ). 

Our helium stars are calculated for a metallicity of Z = 0.02 
(mainly isotopes of C, N, O and Ne). In our default models, we 
assumed a mixing-length para meter of a = l/ H v = 2.0 and 
treated semiconvection as in lLanger et al.1 d 1 9831) adopting an ef¬ 


ficiency parameter qsem = 1 dLanget 1991 ). Convective core 


overshooting was not applied to our models, because the convec¬ 
tive core mass in core helium burning stars increases as function of 
time, which results in a chemical discontinuity at the upper edge of 
the convective core, thereby hindering mixing (cf. Sectio n |3). We 
tested a number of models using a = 1/H p = 1.5 0Langenll99ll) 
which resulted in only slightly smaller final metal-core masses (by 
< 2 per cent). However, a smaller value of a did result in more 
stripping from donor stars in very wide orbits (cf. Section I773l l. In 
our computations, we did not include effects of X-ray irradiation 
of the donor star. Although X-ray irradiation may lead to a cyclic 
mass-transfer behavior (at least for low-mass, hydrogen-rich donor 
stars), previous studies show that these effects are relatively unim¬ 
portant for the final do nor mass and the orbital evolution of X-ray 
binaries with RLO (e.g. lBenvenuto et alj2015l) . 

For evolved helium stars, the typical mass-transfer rates are 
3-4 orders of magnitude larger than the Eddington accretion limit 
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of a NS (~ 4 x 10~ 8 Mq yr -1 , see Section[5]l. This implies that 
> 99.9 per cent of the transferred material is lost from the sys¬ 
tem, presumably in the form of a jet, or a wind from the disk, 
and we assume that it carries the specific orbital angular momen¬ 
tum of the NS, i.e. following the so-called isotropic re-emission 
model fehattacharva & van den I leuvell 199 ll : lTaurisJ 19~%1 ). For re¬ 
cent applications of BEC to X-tay bi naries and further descrip- 


:s ana turther descrtp - 

H fcoiiL I2012LI2013I1 : 


To probe the systems leading to Fe CCSNe, EC SNe and mas¬ 
sive WDs, and to study the effect of mass transfer initiated at differ¬ 
ent evolutionary stages, we investigated systems with helium star 
donors of initial masses of A/H e ,i = 2.5 — 10.0 Mq and initial 
orbital periods, P or b,i = 0.06 — 120 days. For simplicity, we as¬ 
sumed in all cases an initial mass of Mns = 1.35 Mq for the ac¬ 
creting NS. We also evolved isolated helium stars for comparison 
with our binary stars. All configurations are summarized in TableQ] 
(see Section[3]for a description of the parameters). 

The effect of stellar wind mass loss is less important for 
our results. The reason is that orbital dynamics, and more impor¬ 
tantly, the stripping of the donor star, is completely dominated by 
the Roche-lobe overflow. For a compilation of di fferent mass-los s 
rates of Wolf-Rayet (na ked he lium) st ars, se e Yoon et ail d2010h . 
Flere we followed IWe llstein & Langer ( 199f j) an d a pplied the stel¬ 
lar wi nd mass-loss rates of Hamann et al.1 d 19821) or lHamann et al.l 
Jl995h . depending on the luminosities of our relatively low-mass 
helium stars. For the lightest helium stars with an initial mass of 
Mue,i ^ 3.0 Mq we neglected wind mass loss in most of our 
models, except for the widest systems. Applying different prescrip¬ 
tions for the wind mass-loss rate would in effect partly be mimicked 
by choosing a different initial value of P or b- 


3 RESULTS 

In Fig.Q] we show four examples of our modelling with very differ¬ 
ent outcomes. By gradually changing the initial helium star mass, 
Mue,i from 2.6 to 3.5 Mq, and the initial orbital period, P or b,i 
from 0.06 to 2.0 days, the outcome of the resulting remnant varies 
from a CO WD. or an ONeMg WD, to a NS produced via an EC SN 
or an Fe CCSN. In Sections [3. lH3.4l we discuss in detail the results 
of each of these four different models. The final outcome of most 
of our calculations with MHe,i = 2.5 — 3.5 Mq is plotted in Fig.[2] 
In TableQ] we list all relevant parameters for the outcome of 
each calculated system. The parameters listed are the following: 
Mue,i is the initial mass of the naked helium star; P or b,i is the ini¬ 
tial orbital period (always assuming a 1.35 Mq NS companion); 
The RLO Case BA, BB or BC indicates if mass transfer is initi¬ 
ated during core helium burning, helium shell burning, or beyond; 
M C ore,f is the final mass of the metal core of the helium star prior 
to the SN or formation of a WD; is the final mass of he¬ 

lium in the envelope; M*f is the final total mass of the helium star; 
Porb.f is the final orbital period prior to the SN or formation of a 
WD; AMns is the amount of material accreted by the 1.35 Mq 
NS (assuming Eddington-limited accretion, cf. Section [5}: At x is 
the total duration of the X-ray phase, including detached epochs 
in-between; |M^ e ax | is the maximum mass-loss rate of the helium 
star donor during stable Case BA, Case BB or Case BC RLO; |Mf | 
is the mass-transfer rate in our last calculated model. The last three 
columns yield the final fate of the helium star, the merger time of 
the resulting double degenerate system (Section[63} and additional 
comments, respectively. The very first column to the left indicates 


whether or not wind mass loss is included in the evolution of the 
helium star. 


3.1 Evolution leading to an iron core-collapse SN 

Using our stellar evolution code, we were able to evolve the helium 
stars until the onset of silicon burning, a few days prior to the grav¬ 
itational core collapse. In Fig. [3] we have plotted as an example the 
Kippenhahn diagram for the evolution of an initial 3.5 Mq helium 
star with P or b,i = 2.0 days. The mass loss of the helium star is 
initially caused only by its stellar wind (cf. Section[2j. At an age of 
t = 1.426 Myr (at the second hump in the black curve, shortly af¬ 
ter the onset of central carbon burning), the expanding helium star 
fills its Roche lobe with a radius of R = 5.68 Rq and initiates mass 
transfer (Case BC RLO) to the NS at a rate of ~ 10 -4 Mq yr -1 
(see Section |5j. At this stage the remaining stellar lifetime of our 
calculation is log(t* — t) ~ 3.75, i.e. some 5600 yr prior to the 
core collapse. The preceding hump seen in the total mass profile of 
the donor star at log(t* — t) ~ 4.70 is caused by a sudden change 
in the calculated wind mass-loss rate when the luminosity of the 
star increases significantly following the depletion of core helium 
burning. After a sequence of carbon shell burning stages, central 
oxygen burning is ignited about 10 yr prior to core collapse (when 
log(f* — t) ~ 1.0). During the subsequent oxygen shell burning 
stages, silicon is efficiently mixed via convection out to a mass co¬ 
ordinate of ~ 1.55 Mq. This is clearly seen in Fig.[4](upper panel) 
where we have plotted the chemical abundance profile of our final 
stellar model. At this stage (a few days prior to core collapse during 
which silicon burning would take place, producing an iron core) the 
exploding star has a total mass of M*,f = 2.39 Mq with a metal 
core of M core! f = 1.81 Mq. 

In the example illustrated here, there is no doubt that the fi¬ 
nal fate of the donor star is an Fe CCSN. However, in other cases 
it is less trivial to predict the nature of the SN explosion. This is 
the case for stars with a final metal core mass close to the Chan¬ 
drasekhar mass. We could not evolve our stars until the very onset 
of the collapse. However, as argued in Paper I, we can compare 
our evolutionary tracks in the (p c , T c )-plane wit h those availabl e 
in t he literature for iso lated massive stars, e.g. lUmeda et alj 1 201 2l) 
an d[jone s et all j20 1 3l ) . As a simple rule of thumb, we find that evo¬ 
lutionary tracks where the post-carbon burning central temperature, 
T c rises above the value of T c during carbon burning, will eventu¬ 
ally ignite oxygen (typically at T c > 10 9 K, depending on p c ), 
and later burn silicon to produce an iron core and finally undergo 
an Fe CCSN. This is the case for our models which have a metal 
core, M core ,f > 1.43 Mq. We therefore adopt this mass as an 
approximate threshold limit separating Fe CCSNe from EC SNe. 

For low-mass metal cores the exact boundary between Fe CC¬ 
SNe and EC SNe depends on th e locatio n of ignitio n of the off- 
centre neon and oxygen burning dTimmes et al.ll 19941) . and also on 
the propagation of the neon-oxygen flame which is sensitive to 
mixin g in the convective layers across the flame front djones et aD 
120141) . The latter process may affect the electron fraction and the 
density in the central region, and thereby the electron-capture effi¬ 
ciency, and thus the nature of the SN. 


3.2 Evolution leading to an electron-capture SN 

For stars which do not ignite oxygen, but for which their ONeMg 
cores are more massive than 1.37 Mq, the final fate is an EC SN 
dNomotdll987l ; IPodsiadlowski et akll2004l ; iTakahashi et alj|2013l) . 
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Figure 1. Central temperature versus central density for four selected helium star-NS models which terminate their evolution as an Fe CCSN, EC SN, 
ONeMg WD and CO WD, respectively, cf. Sections 13.1113.41 Carbon burning (marked by C) is seen in the first three cases, and oxygen burning (marked by 
O) occurs in the most massive star. Epochs of RLO onset (solid circles) and RLO termination (open squares) are shown along each track. (Note, in some cases 
the first episode of mass transfer occurs for log p c < 4.8.) 
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Figure 2. Sequences of final core mass, M cor e,f as a function of initial helium star mass, Mne,i for four different values of the initial orbital period. Also 
shown (dashed line) are the estimated final core masses for single helium stars. M core f is defined as the metal core (i.e. where the mass fraction of helium is 
less than 10 per cent). The coloured regions indicate the final destiny of the helium stars: light blue (with open triangle symbols) —>■ ONeMg WDs; pink (open 
stars) —»• EC SNe; and purple (solid stars) —> Fe CCSNe. In all cases the accreting companion star is (initially) a 1.35 Mq NS. 
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Table 1. Stellar and binary parameters of 68 helium star-NS systems, and 8 single helium stars, evolved to the pre-SN (or WD) stage - see Section[3] 


Wind 

M H e,i 

(%) 

-Porb.i 

(days) 

RLO 

Case 

M core,f 

(Mq) 

A/f env 

M He,f 

(Mq) 

f 

(Mq) 

P orb,f 

(days) 

am ns 

(Mq) 

A t x 

(yr) 

|Mf“| 

(M 0 yr -1 ) 

|M f | 

(M 0 yr -1 ) 

final fate 

Tgrw 

(Myr) 

comment 

yes 

2.5 

_ 

_ 

1.29 

1.099 

2.41 

_ 

- 

_ 

_ 

_ 

ONeMg 

_ 

single star 

yes 

2.6 

- 

- 

1.37 

0.961 

2.37 

- 

~ 

- 

- 

- 

EC-SN 

- 

single star 

yes 

2.7 

- 

- 

1.41 

1.010 

2.46 

- 

- 

- 

- 

- 

EC-SN 

- 

single star 

yes 

2.8 

- 

- 

1.46 

0.984 

2.49 

- 

- 

- 

- 

- 

FeCCSN 

- 

single star 

yes 

2.9 

- 

- 

1.50 

0.982 

2.54 

- 

- 

- 

- 

- 

FeCCSN 

- 

single star 

yes 

3.0 

- 

- 

1.58 

0.970 

2.60 

- 

- 

- 

- 

- 

FeCCSN 

- 

single star 

yes 

3.2 

- 

- 

1.70 

0.973 

2.72 

- 

- 

- 

- 

- 

FeCCSN 

- 

single star 

yes 

3.5 

- 

- 

1.83 

1.003 

2.91 

- 

- 

- 

- 

- 

FeCCSN 

- 

single star 

yes 

3.0 

120 

- 

1.58 

0.968 

2.59 

145.5 

- 

- 

- 

- 

FeCCSN 

l.ell 

no RLO at all 

yes 

3.0 

100 

BC 

1.57 

0.750 

2.38 

109.7 

4.8e-5 

1.2e3 

2.7e-4 

2.7e-4 

FeCCSN 

5.el0 


yes 

3.0 

80 

BC 

1.57 

0.736 

2.37 

87.1 

5.2e-5 

1.3e3 

2.6e-4 

1.4e-4 

FeCCSN 

3.el0 


yes 

3.0 

50 

BC 

1.58 

0.689 

2.32 

53.2 

6.3e-5 

1.5e3 

2.7e-4 

l.le-3 

FeCCSN 

6.6e9 


yes 

2.5 

20 

BC 

1.29 

0.137 

1.45 

20.3 

6.1e-4 

1.4e4 

2.2e-4 

6.7e-5 

ONeMg 

l.le8 


yes 

2.6 

20 

BC 

1.37 

0.189 

1.56 

19.3 

6.3e-4 

1.5e4 

8.8e-5 

8.1e-5 

EC-SN 

1.9e8 


yes 

2.7 

20 

BC 

1.42 

0.150 

1.59 

19.7 

4.5e-4 

l.le4 

2.2e-4 

7.7e-5 

EC-SN 

1.9e8 


yes 

2.8 

20 

BC 

1.47 

0.217 

1.72 

19.4 

3.2e-4 

7.8e3 

3.0e-4 

detached 

FeCCSN 

2.1e8 


yes 

2.9 

20 

BC 

1.52 

0.342 

1.91 

19.3 

2.2e-4 

5.4e3 

3.2e-4 

3.9e-8 

FeCCSN 

2.6e8 


yes 

3.0 

20 

BC 

1.58 

0.548 

2.17 

20.0 

1.4e-4 

3.5e3 

3.0e-4 

8.1e-5 

FeCCSN 

3.8e8 


yes 

3.2 

20 

BC 

1.70 

0.919 

2.67 

23.8 

2.8e-5 

7.0e2 

1.6e-4 

1.6e-4 

FeCCSN 

1.0e9 


yes 

3.5 

20 

- 

1.83 

1.007 

2.91 

25.7 

- 

- 

- 

- 

FeCCSN 

1.4e9 

no RLO at all 

yes 

3.0 

10 

BC 

1.57 

0.375 

2.00 

9.51 

2.3e-4 

5.6e3 

2.9e-4 

4.e-10 

FeCCSN 

4.1e7 


yes 

3.0 

5.0 

BC 

1.57 

0.246 

1.86 

4.61 

4.1e-4 

1.0e4 

2.0e-4 

detached 

FeCCSN 

4.7e6 


no 

2.5 

2.0 

BC 

1.29 

0.093 

1.40 

1.94 

6.8e-4 

2.2e4 

2.5e-4 

1.5e-5 

ONeMg 

2.2e5 


no 

2.6 

2.0 

BC 

1.37 

0.079 

1.46 

1.78 

8.3e-4 

2.2e4 

2.9e-4 

3.8e-5 

EC-SN 

2.8e5 


no 

2.7 

2.0 

BC 

1.42 

0.092 

1.53 

1.62 

6.9e-4 

1.8e4 

2.5e-4 

3.7e-6 

EC-SN 

2.2e5 


no 

2.8 

2.0 

BC 

1.47 

0.121 

1.62 

1.47 

6.2e-4 

1.6e4 

2.0e-4 

4.2e-5 

FeCCSN 

1.8e5 


no 

2.9 

2.0 

BC 

1.52 

0.152 

1.72 

1.34 

5.5e-4 

1.4e4 

1.7e-4 

3.9e-5 

FeCCSN 

1.5e5 


no 

3.0 

2.0 

BC 

1.57 

0.176 

1.80 

1.23 

4.8e-4 

1.2e4 

1.8e-4 

2.8e-6 

FeCCSN 

1.2e5 


yes 

3.2 

2.0 

BC 

1.66 

0.244 

1.98 

1.71 

4.0e-4 

9.7e3 

1.5e-4 

7.7e-8 

FeCCSN 

3.4e5 


yes 

3.5 

2.0 

BC 

1.81 

0.488 

2.39 

1.78 

2.6e-4 

6.5e3 

1.9e-4 

l.le-7 

FeCCSN 

5.4e5 

Section|3.1 1 

yes 

3.7 

2.0 

BC 

1.93 

0.754 

2.75 

2.10 

1.8e-4 

4.5e3 

1.6e-4 

1.5e-4 

FeCCSN 

l.le6 


yes 

4.0 

2.0 

BC 

2.02 

0.850 

2.95 

2.45 

1.5e-4 

3.6e3 

1.0e-4 

8.9e-5 

FeCCSN 

1.9e6 


yes 

5.0 

2.0 

- 

2.19 

0.861 

3.15 

3.91 

- 

- 

- 

- 

FeCCSN 

6.4e6 

no RLO at all 

no 

3.2 

2.0 

BC 

1.70 

0.232 

1.98 

1.05 

3.5e-4 

8.7e3 

2.2e-4 

4.7e-5 

FeCCSN 

8.5e4 


no 

3.5 

2.0 

BC 

1.91 

1.00 

2.99 

1.30 

1.8e-4 

4.4e3 

2.8e-4 

1.6e-3 

FeCCSN 

4.5e5 


no 

2.5 

0.5 

BB 

1.27 

0.053 

1.41 

0.500 

1.0e-3 

4.2e4 

1.0e-4 

2.3e-5 

ONeMg 

5960 


no 

2.6 

0.5 

BB 

1.35 

0.003 

1.35 

0.475 

9.0e-4 

3.9e4 

1.2e-4 

2.8e-4 

ONeMg 

4940 


no 

2.7 

0.5 

BB 

1.41 

0.048 

1.48 

0.415 

8.4e-4 

3.2e4 

1.4e-4 

7.7e-5 

EC-SN 

5430 

Section|3.2| 

no 

2.8 

0.5 

BB 

1.46 

0.072 

1.56 

0.376 

8.2e-4 

2.9e4 

1.6e-4 

8.1e-4 

FeCCSN 

4350 


no 

2.9 

0.5 

BB 

1.51 

0.090 

1.64 

0.341 

7.8e-4 

2.6e4 

1.8e-4 

1.4e-7 

FeCCSN 

3480 


no 

3.0 

0.5 

BB 

1.56 

0.112 

1.73 

0.311 

7.2e-4 

2.2e4 

1.8e-4 

6.6e-5 

FeCCSN 

2880 


yes 

3.2 

0.5 

BB 

1.66 

0.142 

1.86 

0.387 

6.1e-4 

1.7e4 

1.7e-4 

detached 

FeCCSN 

5270 


yes 

3.5 

0.5 

BB 

1.80 

0.183 

2.07 

0.365 

5.0e-4 

1.2e4 

2.0e-4 

2.4e-5 

FeCCSN 

4860 


yes 

10.0 

0.5 

BC 

2.17 

0.446 

3.02 

2.94 

2.7e-5 

9.3e2 

5.7e-5 

2.4e-8 

FeCCSN 

2.6e6 


no 

3.2 

0.5 

BB 

1.68 

0.134 

1.86 

0.261 

6.3e-4 

1.8e4 

2.7e-4 

9.3e-6 

FeCCSN 

1770 


no 

3.5 

0.5 

BB 

1.88 

0.203 

2.15 

0.207 

4.9e-4 

1.2e4 

3.9e-4 

1.0e-4 

FeCCSN 

1050 


no 

2.5 

0.1 

BB 

1.22 

0.015 

1.24 

0.108 

4.2e-3 

1.2e5 

4.1e-5 

2.1e-5 

ONeMg 

104 


no 

2.6 

0.1 

BB 

1.28 

0.004 

1.29 

0.098 

3.6e-3 

l.le5 

5.3e-5 

l.le-4 

ONeMg 

76.8 

Section|3.3| 

no 

2.7 

0.1 

BB 

1.34 

0.004 

1.35 

0.088 

3.1e-3 

9.2e4 

6.6e-5 

7.4e-4 

ONeMg 

55.5 


no 

2.8 

0.1 

BB 

1.39 

0.018 

1.43 

0.079 

2.7e-3 

8.2e4 

8.4e-5 

8.4e-3 

EC-SN 

61.9 


no 

2.9 

0.1 

BB 

1.45 

0.033 

1.50 

0.070 

2.1e-3 

7.3e4 

l.le-4 

l.le-2 

FeCCSN 

44.8 

Paper I 

no 

3.0 

0.1 

BB 

1.50 

0.040 

1.58 

0.063 

2.0e-3 

6.5e4 

1.3e-4 

1.9e-2 

FeCCSN 

35.3 


yes 

3.2 

0.1 

BB 

1.60 

0.054 

1.71 

0.056 

1.7e-3 

5.3e4 

1.7e-4 

4.6e-4 

FeCCSN 

26.5 


yes 

3.5 

0.1 

BB 

1.76 

0.077 

1.88 

0.048 

1.2e-3 

4.1e4 

2.2e-4 

detached 

FeCCSN 

17.2 


yes 

3.7 

0.1 

BB 

1.86 

0.079 

2.00 

0.047 

l.le-3 

3.4e4 

2.2e-4 

detached 

FeCCSN 

16.4 


yes 

4.0 

0.1 

BB 

1.96 

0.094 

2.12 

0.048 

9.4e-4 

2.9e4 

2.1e-4 

detached 

FeCCSN 

17.6 


yes 

5.0 

0.1 

BB 

2.10 

0.112 

2.36 

0.062 

6.6e-4 

1.9e4 

l.le-4 

6.0e-3 

FeCCSN 

39.5 


yes 

6.0 

0.1 

BB 

2.15 

0.109 

2.37 

0.064 

6.7e-4 

1.9e4 

1.0e-4 

3.0e-5 

FeCCSN 

40.2 


no 

3.2 

0.1 

BB 

1.61 

0.057 

1.71 

0.052 

1.7e-3 

5.4e4 

2.0e-4 

4.6e-4 

FeCCSN 

21.3 


no 

3.5 

0.1 

BB 

1.81 

0.074 

1.93 

0.039 

1.0e-3 

3.9e4 

3.9e-4 

detached 

FeCCSN 

9.76 


no 

2.5 

0.08 

BB 

1.21 

0.017 

1.24 

0.085 

4.6e-3 

1.3e5 

4.7e-5 

detached 

ONeMg 

55.1 


no 

2.6 

0.08 

BB 

1.25 

0.097 

1.27 

0.078 

4.2e-3 

1.2e5 

5.9e-5 

3.4e-5 

ONeMg 

42.6 


no 

2.7 

0.08 

BB 

1.32 

0.005 

1.32 

0.071 

3.6e-3 

l.le5 

7.5e-5 

1.9e-4 

ONeMg 

31.7 


no 

2.8 

0.08 

BB 

1.38 

0.005 

1.38 

0.064 

2.9e-3 

9.1e4 

9.6e-5 

2.8e-3 

EC-SN 

32.7 


no 

2.9 

0.08 

BB 

1.41 

0.016 

1.45 

0.057 

2.6e-3 

8.5e4 

1.2e-4 

3.2e-4 

EC-SN 

25.7 


no 

3.0 

0.08 

BB 

1.45 

0.021 

1.49 

0.052 

2.3e-3 

8.0e4 

1.6e-4 

6.0e-2 

FeCCSN 

19.9 


yes 

3.2 

0.08 

BB 

1.58 

0.046 

1.68 

0.043 

1.6e-3 

5.9e4 

2.3e-4 

2.3e-3 

FeCCSN 

12.9 


yes 

3.5 

0.08 

BB 

1.71 

0.055 

1.81 

0.035 

8.5e-4 

5.1e4 

4.3e-4 

detached 

FeCCSN 

7.23 


no 

2.5 

0.06 

BA 

0.80 

0.036 

0.87 

0.102 

3.3e-2 

5.1e7 

5.9e-5 

detached 

CO 

127 


no 

2.6 

0.06 

BA 

0.76 

0.045 

0.84 

0.104 

2.5e-2 

4.9e7 

7.9e-5 

detached 

CO 

140 

Section|3.4| 

no 

2.7 

0.06 

BA 

0.67 

0.064 

0.78 

0.118 

1.3e-3 

6.0e4 

l.le-4 

detached 

CO 

220 


no 

2.8 

0.06 

BA 

0.0 

1.79 

2.22 

? 

l.le-3 

4.7e4 

1.6e-4 

2.3e-4 

CE(?) 

- 

fate? 

no 

2.9 

0.06 

BA 

0.0 

1.88 

2.30 

? 

9.3e-4 

4.3e4 

2.3e-4 

9.5e-4 

CE(?) 

- 

fate? 

no 

3.0 

0.06 

BA 

0.0 

2.13 

2.30 

? 

8.0e-4 

3.9e4 

2.6e-4 

4.4e-4 

CE(?) 

- 

fate? 

yes 

3.2 

0.06 

BA 

0.0 

- 

2.67 

- 

6.7e-4 

3.4e4 

3.6e-4 

988. 

CE 

- 

merger 

yes 

3.5 

0.06 

BA 

0.0 

- 

3.06 

- 

5.2e-4 

2.8e4 

l.le-3 

1004. 

CE 

- 

merger 


© 0000 RAS, MNRAS 000, 000-000 






6 Tauris, Langer & Podsiadlowski 



log (t.-t)/yr 



Gfl 

O 


Figure 3. Kippenhahn diagram of the 3.5 M© helium star plotted in Fig.0 
and leading to the final structure shown in the upper panel of Fig. \4\ The 
plot shows cross-sections of the helium star in mass coordinates from the 
center to the surface of the star, along the y-axis, as a function of stellar age 
on the rr-axis. The value {t * — t) is the remaining time of our calculations, 
spanning a total time of f* = 1.433 Myr. The green hatched areas denote 
zones with convection; red colour indicates semiconvection. The intensity 
of the blue/purple colour indicates the net energy-production rate. Carbon 
burning was ignited at log(t* — t) « 4.0 and oxygen burning was ignited 
at log(t* — t) « 1.0. Our calculations ended just before silicon ignition, a 
few days prior to the gravitational collapse. 


We find that ten of our models (including six with P or b ^ 
2.0 days, cf. Fig. [2} produce EC SNe. All these stars have metal 
cores with a mass, M cor e,f — 1-37 — 1.43 M©. 

In the central panel of Fig. [4] we have plotted the chemical 
abundance profile of our final stellar model resulting from a helium 
star donor with Mne,i = 2.7 M© and P or b,i = 0.5 days. The plot 
shows the chemical composition shortly before the loss of pressure 
support, due to inverse /3-decay, leading to the EC SN event. In this 
particular example, the exploding star has an envelope containing 
less than 5 x 10 -2 M© of helium. Therefore, this EC SN would 
most likely be classified as a Type Ic SN, see discussion in Sec¬ 
tion [472] 


3.3 Evolution leading to an ONeMg WD 

In the lower panel of Fig. [4] we have plotted the chemical abun¬ 
dance profile of a stripped donor star whose progenitor was a 
2.6 M© helium star with P or b,i = 0.1 days. Following Case BB 
RLO, the final mass is 1.29 M© and this star will leave behind an 
ONeMg WD. It is interesting to note the pureness of this ONeMg 
WD. Almost the entire WD consists of only three elements: oxygen 
(38.0 per cent), neon (48.9 per cent) and magnesium (10.2 per cent). 
It is surrounded with a tiny envelope of less than 4 x 10 -3 M© he¬ 
lium and < 3 x 10 ~ 2 M© carbon . If such a WD were to undergo 
pulsations dWinget & Kepleill20081) , one might be able to verify its 
structure observationally. 


3.4 Evolution leading to a CO WD 

Last but not least, we present the evolution of a similar 2.6 M© he¬ 
lium star, but now placed in a tighter orbit with P or b,i = 0.06 days. 
This star fills its Roche lobe already during early core helium burn¬ 
ing. The resulting mass transfer causes the star to lose its outer 
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Figure 4. Final chemical abundance structure of three evolved helium stars. 
Upper panel: a helium star with an initial mass of 3.5 M© placed in a 
binary with a 1.35 M© NS and an orbital period of 2.0 days. The final mass 
of the exploding star is 2.39 M©, with a metal core of 1.81 M©, leading 
to an Fe CCSN of Type lb. Our last calculated model #2680 is just before 
the onset of silicon burning, a few days prior to the explosion, at a helium 
star age of t = 1.433 Myr. Central panel: a helium star with an initial 
mass of 2.7 M© placed in a binary with a 1.35 M© NS and an orbital 
period of 0.5 days. After Case BB RLO the final mass of the exploding star 
is 1.48 Mq, with a metal core of 1.41 Mq, leading to an EC SN, probably 
of Type Ic (see text). Three magnesium isotopes ( 24 Mg, 25 Mg, 26 Mg) are 
plotted as a solid, dashed and a dotted green line, respectively, yielding a 
total mass fraction of 10.9 per cent. Lower panel: a helium star with an 
initial mass of 2.6 M© placed in a binary with a 1.35 M© NS and an 
orbital period of 0.1 days. The final remnant is an almost pure ONeMg WD 
with a mass of 1.29 M©. The total mass fraction of Mg is 10.2 per cent. 


© 0000 RAS, MNRAS 000, 000-000 
















































Ultra-stripped SNe: progenitors and fate 


7 



Figure 5. Kippenhahn diagram of models 1 to 554 of the 2.6 Mq helium 
star plotted in Figs. and [6] The energy production is indicated in blue 
shades. The hatched green colour indicates convection. The central helium 
burning is completely shut-off as a result of the intense mass loss from the 
outer layers of the donor star during the first mass-transfer phase (Case BA 
RLO). Central helium burning is only reignited after the RLO has termi¬ 
nated and the central temperature has risen again. See also Fig. [7] 



Figure 7. Central temperature versus central mass density for the donor 
star plotted in Figs. |T| [%] and [6] leading to a 0.84 Mq CO WD. The red 
lines indicate rates of energy production from central helium burning (the 
difference between adjacent lines is one order of magnitude). Shortly into 
the first mass-transfer episode (Case BA), central helium burning is tem¬ 
porarily shut off as a result of the decreasing core temperature. The reason 
for this is work required to expand the inner region while replacing the outer 
layers of the star which are lost via RLO. 


~ 1.6 Mq of material in less than 50 000 yr, cf. Fig. [5] As a di¬ 
rect consequence of this mass loss, the star is driven out of thermal 
equilibrium. In order for the star to replace the lost envelope with 
material from further below and to remain in hydrostatic equilib¬ 
rium, an endothermic expansion of its inner region (requiring work 
against gravity) causes the surface luminosity to decrease by more 
than 3 orders of magnitude, cf. Fig. [ 6 ] This effect has an important 
consequence for the nuclear burning in the star. The core expan¬ 
sion causes the central temperature to decrease (Fig.0, and given 
the extreme dependence of the triple-a process on temperature, the 
result is that all fusion processes are completely quenched for about 
80 000 yr. Thereafter, when the central temperature has risen again, 
core helium burning is reignited. 



Age (Myr) 

Figure 8. Mass-transfer rate, \l\d 2 \ (blue peaks) and orbital period, P Q rb 
(green curve) as a function of stellar age for the full evolution of the 
2.60 Mq helium star orbiting a NS with an initial orbital period of 
0.06 days (also shown in Figs. 00 There are three phases of mass trans¬ 
fer: Case BA (onset during core helium burning); Case BAB (a subsequent 
episode of mass transfer when the remaining part of the donor star expands 
during helium shell burning); and finally mass transfer as an UCXB when 
the detached CO WD fills its Roche lobe as a result of strong loss of orbital 
angular momentum via gravitational wave radiation. The dashed red line 
indicates the Eddington accretion rate of the NS (~4 x 10 —8 Mq yr -1 
for accretion of helium). The evolution of the orbital period is dominated 
by either mass transfer (RLO) or gravitational wave radiation. 


In Fig. [ 8 ] we have plotted the full evolution of the 2.6 Mq 
helium star (with P OI b,i = 0.06 days) as a function of age. There 
are three phases of mass transfer: Case BA was described in the 
previous section; Case BAB is a result of stellar expansion during 
helium shell burning; and finally the star fills its Roche lobe again, 
130 Myr after it has become a 0.84 M 0 CO WD, as a result of 
orbital shrinkage due to gravitational wave radiation. At this point 
P orb = 32 sec and the system becomes an ultra-compact X-ray 
binary (UCXB). However, as a result of the relatively large mass 
ratio between the WD and the NS, this UCXB is short live d and 
most likely leads to a merger event dvan Haaften et al.ll2012h . De¬ 
pending on the NS equation-of-state, the outcome of this merger 
is either a massive NS or a black hole ( possi bl y acco mpanied with 
a 7 -ray burst-like event, see e.g. iDermer & Atovanll2006h . It has 
been shown that a significant fraction (50 to 80 per cent ) of the 
WD mass could be ejected during the merger process llMetzgerl 
120121 ) ; however t he details of the mode lling of such events still 
remain uncertain IPasch alidis et af1 l201lli . Even if assuming fully 
conservative mass transfer, the gravitational mass of the resulting 
compact object formed by the merger considered here is at most 
~ 2.1 M 0 , i.e. close to the value of the recently measu red NS 
mass of ~ 2.01 ± 0.04 Mq for PSR J0348+0432 dAntoniadis et al.l 
l2013lf . Therefore, we conclude that in this case the final product is 
most likely still a NS. 


3.5 Evolution of massive helium star-NS binaries 

We explored the evolution of a 10 M 0 helium star orbiting a NS 
with an initial orbital period of either P or b,i = 0.1, 0.5 or 2.0 days. 
Cygnus X-3 is a potential observational example of such a system. 
It contains a compact obj ect with a Wolf-Rayet star companion 
dvan Kerkwiik et al.ll 19921) and has P OIb = 4.8 hours. The exact 
nature_of the com pact object (NS or BH), however, is still unknown 
dZdziarski et al.ltloTll . and references therein). 


© 0000 RAS, MNRAS 000, 000-000 








































8 Tauris, Langer & Podsiadlowski 



!°g T eff (K) 

Figure 6. HR-diagram of the 2.6 Mq helium star donor discussed in Section 13741 The evolutionary calculations start on the helium zero-age main sequence 
(ZAMS, model 1 marked with a black square). Three phases of RLO mass transfer are indicated with blue colour along the evolutionary track. Central helium 
burning is temporarily shut off during the first RLO. It is reignited shortly after model 537, and at model 694 (f = 6.8 Myr) central helium is depleted, 
yielding a CO core made of 42 per cent carbon and 56 per cent oxygen. The subsequent expansion of the star during helium shell burning causes a second 
phase of mass transfer. Thereafter, the 0.84 Mq CO WD settles on the cooling track before gravitational wave radiation causes it to fill its Roche lobe again 
(t — 130 Myr) at P or b = 32 sec, leading to a short UCXB-phase before the system merges (model 1938, see text). The inset is a zoom-in of the area where 
the track crosses itself near the location of model 1 (f = 0). 


The orbital evolution, and thus the fate, of such a system is 
determined by the competition between gravitational wave radia¬ 
tion and stellar wind mass loss from the massive helium (Wolf- 
Rayet) star, which causes the orbit to shrink and widen, respec¬ 
tively. A lso spin-o rbit cou pling may become important in such 
systems foetmers et al.ll2008l) . which is, however, not considered 
here. For P OI b,i = 0.1 days, the orbital damping by gravita¬ 
tional wave radiation is efficient and causes the helium star to ini¬ 
tiate RLO at P or b = 0.073 days while it is still rather massive 
(about 8.90 Mq). As a result of the high mass ratio between the 
helium star donor and the accreting NS at this point, the system 
becomes dynamically unstable, enters a CE phase, and merges. For 
Porb,i = 2.0 days, and with our applied wind mass-loss prescrip¬ 
tion, the system widens significantly and reaches a large dimen¬ 
sion, preventing the radius of the helium star from filling its Roche 
lobe at any time. For orbital periods in-between, there are solutions 
which result in stable RLO. We find, for example, a solution for 
Porb,i = 0.5 days (Table[]J. 

For all such massive helium star systems which undergo stable 
RLO, the mass-transfer phase is very short lasting (< 10 li yr) since 
these stars are evolved to an advanced stage, close to carbon shell 
burning, before they expand sufficiently to compensate for the or¬ 
bital widening by wind mass loss and initiate RLO. Hence, their re¬ 
maining nuclear lifetime is short before they collapse and explode 
(a situation which is similar for low-mass helium stars initiating 
RLO in very wide-orbit systems, cf. Section EOl i. For the particular 
case with P or b i = 0.5 days, our adopted stellar wind mass loss 
jHamann et alji 19951 ) reduces the mass of the helium star down to 


3.04 Mq prior to RLO. During the subsequent mass transfer, the 
evolved helium star only loses 0.02 Mq before it explodes. 

We caution that different applied wind mass-loss rates of mas¬ 
sive Wolf-Rayet stars can result in significantly different pre-SN 
stellar masses . For example, reducing the wind mass-loss rate of 
[Hamann^tal. i 1V95) by a factor of 10, or applying the prescription 
of Nugis & Lamersl ' 2000h . yields a pre-SN mass of about 7.0 Mq 
in the above-mentioned case. Therefore, we notice that our helium 
stars models with Mu e ,i = 5.0, 6.0 and 10.0 Mq, which were 
found to have very similar metal core masses between 2.10 and 
2.20 Mq by the time of their explosion, could be more massive 
if the wind mass-loss rate is reduced. For a discussion and a com- 
parison between various prescriptions of wind mass loss-rates, see 
Yoon et alj d20 1 dl) . 


4 OBSERVATIONAL PROPERTIES OF 
ULTRA-STRIPPED SUPERNOVAE 

4.1 Remaining helium envelope mass prior to the SN 

The modelling of synthetic SN spectra by lHachinger et alj d2012h 
suggests that more than 0.06 Mq of helium is needed for helium 
lines to become visible in optical/IR spectra, and thus for classifica¬ 
tion as a Type lb rather than a Type Ic SN event. However, this ques¬ 
tion also depends on the amount of nickel synthesized during the 
explosion and how w ell it is m ixed into t he helium-rich material of 
the ejected envelope foessart et al.ll201 lh . given that the detectable 
helium is excited by the radioactive decay of nickel. Furthermore, 
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Initial He—star mass, M Hei (M 0 ) 

Figure 9. Total amount of helium prior to explosion for the models from 
Fig. [2] leading to EC SNe (open stars) and Fe CCSNe (solid stars). The 
solid lines connect sy stems with equal values of P 0 rb,i- According to 
lHachinger et dl 120131 . ejection of at least 0.06 Mq of helium is needed 
to detect a SN as Type lb (yellow and green region). Below this limit the 
SN is likely to be classified as a Type Ic (light blue region). However, this 
also depends on the amount of nickel synthesized - see text. 

explosive helium burning during the S N explosion may also play a 
role in reducing the amount of helium dWooslev & Weaveril 19951) . 

In Fig. [9] we show the amount of helium, in the en¬ 

velopes of most of our pre-SN models with Mu e ,i ^ 3.5 Mq. We 
find that, in general, EC SNe in binaries (particularly close bina¬ 
ries) are most likely to be observed as ultra-stripped Type Ic SNe, 
whereas our resulting Fe CCSNe can be observed as both ultra- 
stripped (faint) Type Ic and Type lb SNe, as well as more regu¬ 
lar Type lb SNe (the latter still with relatively small helium ejecta 
masses, typically less than 0.5 Mq). 

In comparison to the remaining amount of helium, M^f, we 
find from our binary models that the total envelope mass M enVj f = 
M*f — M core , f (i.e. all material outside the CO core) is larger by 
0.03 Mq on average. This additional material is mainly composed 
of carbon which was mixed into the envelope. 


4.1.1 Definition of ultra-stripped SNe 

In Fig. [TO] we show the total amount of envelope mass, M enVj f for 
most of our pre-SN models as a function of the final core mass 
of the exploding star, M core ,f. The dotted line above M env ,f ~ 
0.2 Mq shows the minimum envelope mass ex pected fo r the stars 
which explode first in close massive binaries I lYoon et al.ll2010l) . 
These stars lose mass to a main-sequence (or slightly evolved) com¬ 
panion star via stable RLO prior to the SN. 

In the following we assume that the envelope mass surround¬ 
ing the pre-SN metal core can be taken as a rough estimate of the 
amount of material ejected during the SN explosion. (The deter¬ 
mination of the exact location of the mass cut requires detailed 
computations of the explosion event.) We therefore suggest to de¬ 
fine ultra-stripped SNe as exploding stars whose progenitors are 
stripped more than what is possible with a non-degenerate compan¬ 
ion. In other words, ultra-stripped SNe are exploding stars which 
contain envelope masses < 0.2 Mq and having a compact star 
companion. The compact nature of their companion allows for ex¬ 
treme stripping in a tight binary, in some cases yielding an envelope 
mass ;C 0.008 Mq. 



1.3 1.4 1.5 1.6 1.7 1.8 1.9 


Final core mass, M coref (M 0 ) 

Figure 10. Total envelope mass versus final core mass prior to the explosion 
for the models shown in Figs.[2]an d|9] EC SNe (limited to the grey-shaded 
area) and Fe CCSNe are plotted with open and solid stars, respectively. 
The dotted line shows the remaining envelope mass in pre-SN stars which 
explode first in massive binaries (i.e. after RLO to a main-sequence star 
companion) according to the calculations bv lYoon et all I20K1 ). Hence, we 
use this division to define ultra-stripped SNe as exploding stars in compact 
binaries which have M env f < 0.2 Mq - see text. 

In addition to forming low-mass (1.1 —1.4 Mq) NSs from (al¬ 
most) naked metal cores slightly above the Chandrasekhar niastQ 
we also expect ultra-stripped SNe to produce more massive NSs. 
From Table[T] it is seen that ultra-stripped SNe occur in metal cores 
as massive as 2.15 Mq (cf. the model with Mne.i = 6.0 Mq and 
Torb.i = 0.1 days, and which leads to a post-Case BB RLO stellar 
mass of M*,f = 2.37 Mq and a helium and carbon-rich envelope 
of only M env ,f = 0.22 Mq). Therefore, it seems plausible that 
ultra-stripped SNe may produce NSs with a wide range of masses 
between 1.10 — 1.80 Mq, see Section [6!2l for further discussion. 

We emphasize that ultra-stripped SNe only constitute a small 
fraction of all SNe Ic (< 1 per cent, Paper I), and also note that 
the ejected envelopes of exploding massive Wolf-Rayet stars may 
contain considerable amounts of carbon and oxygen. Whereas the 
pre-SN orbital periods of ultra- stripped SN e arejypically between 
1 hr and 2 days, the models of lYoon et ali d201C)l~) resulted in pre- 
SN orbital periods, at the moment of the first explosion in close 
massive binaries experiencing stable RLO, of at least ~ 10 days. 

As also seen from Figs.l9landllOI besides core mass, the final 
amount of helium and the total envelope mass is also correlated 
with the initial orbital period of the binary prior to Case BB RLO. 
Therefore, we also find that and M env ,f are correlated with 

the final orbital period, P or b,f just prior to the SN, cf. TableQ] 

4.2 Expected light curve properties 

The relatively recent advent of all-sky, high-cadence transient sur¬ 
veys has led to the discovery of many rapidly evolving, luminous 
transients (see iDrout et al .1 uoTil . for details). iDrout et ahl d2014l) 
have estimated the rate of rapid transients (whose luminosity de¬ 
creases by at least 50 per cent in 12 days) to be 4-8 per cent of the 
core-collapse rate (at a redshift z = 0.2). Observationally, these 
newly discovered luminous transients form a heterogeneous group 
of objects, possibly suggesting a variety of explosion mechanisms. 

1 Timm es et a j] h996l) : IWooslev et alj J2002l) . 
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For comparison, the rate of less lumino us SN 2005 ek-lik e 
events is at least 1-3 per cent of the SN la rate jDrout et~aill20 1 3l) . 
As we have argued in Paper I, the oxygen-rich. Type Ic supernova 
SN 2005ek is probably consistent with an ultra-stripped SN, i.e. a 
helium star stripped of most of its heli um by binary interactions, as 
predicted by our models (see also lDrout et al.l2013i) . In Paper I, we 
estimated that the fraction of ultra-stripped SNe to all SNe could be 
as high as 0.1-1 per cent. 

It is beyond the scope of this paper to compute the SN ex¬ 
plosions from our progenitor models. However, the corresponding 
events can be expected to be fast and, in most cases, probably quite 
faint. Based on the photon diffusion time through a homologously 
expan ding SN envelope dArnetd 1 1979L 11982 : iKleiser & Kasenl 
boil . we estimate the rise time r r from the light diffusion time 
through the SN as: 

r r = 5.0 days M^ 4 k 0 ( 2 E ^ 0 1/4 , (1) 

which for the parameter range of our ultra-stripped SN progenitor 
models results in rise times between 12 hr and 8 days (cf. Table[2]l. 
Here, Mo.i is the ejecta mass in units of 0.1 Mg, E$o the SN 
kinetic energy in units of 10 50 erg, and ko.i the opacity in units 
of 0.1 cm 2 g -1 . Here we assumed reo.i = 1, which is suitable for 
electron scattering in singly ionized helium (although the opacity of 
the SN ejecta does depe nd on the temperature a nd the composition 
of the expelled material. lKleiser & Kasenlboi4l . 

The peak b rightness can be estimated according to Arnett’s 
rule dArnetdl 1982h . which results in absolute magnitudes of M = 
— 11.8, —13.4, —15.0, and —16.6 for nickel masses of 0.001, 
0.005, 0.01, and 0.05 Mg, respectively, where we gauged the rela¬ 
tion with the models of Paper I. As the diffusion time is shorter than 
the radioactive decay timescale, the light curves should follow the 
radioactive decay law immediately after the light curve peak, un¬ 
til the SN becomes optically thin to gamma-rays at time ra, when 
the decay of the optical light curve should accelerate. By assuming 
a constant density in the SN ejecta, and a constant velocity of the 
ejecta of v exp = (2EsN/A3ej)°'°, we can estimate the time when 
the electron scattering optical depth of the SN light curve becomes 
equal to one as: 

ra = 25 days M 0 .i Kg /2 E^ 0 1/2 , (2) 


which results in timescales of the order of 1 to 50 days for our 
models (cf. Tabled 

Fast and faint SNe have been suggested to occur also 
based on other models: t hermo nuclear explosions jBildsten et al.l 
120071: iDessart & Hilfiei] l2015il . accretion-induced collapse of 
mass ive WDs j jPessart et al.1 20061 ) as well as core collapse 
dKleiser & Kaseii 2014l t. where for some of them it remains to be 
seen whether they are produced by Nature. As the observed short 
orbital period NS systems give evidence to the scenario elaborated 
here (cf. Section[6]l, we hope that the spectroscopic features of these 
ultra-stripped SNe — which remain to be worked out — will al¬ 
low an unambiguous identification of these events which could then 
serve to measure the NS merger rate directly (Section[6j5]l. 


4.3 Location of ultra-stripped SNe in host galaxies 

To further help identifying optical transients as ultra-stripped SNe 
it is important to know their expected location with respect to star¬ 
forming regions in host galaxies. An ultra-stripped SN is the second 
SN in a binary and thus the travel distance of the system, d, with 
respect to its original birth location as a ZAMS binary, is the prod¬ 
uct of the systemic velocity resulting from the first SN, v sys , and 


Table 2. Expected properties of the light curves resulting from ultra- 
stripped SNe. The assumed explosion energy (2?sn) is given in the left 
column. For each of the different SN ejecta masses (M e j) in the following 
four columns, estimated values of the rise time (r r ) and the decay time (ra) 
of the SN light curve are stated. 


-Esn (erg) 

M e j — 0.2 Mq 

0.1 Mq 

0.03 Mq 

0.01 Mq 

10 50 

r r = 8.4 days 
ra = 50 days 

5.0 days 
25 days 

2.0 days 
7.5 days 

0.9 days 

2.5 days 

3 x 10 50 

T r = 6.4 days 

Td = 29 days 

3.8 days 

14 days 

1.5 days 
4.3 days 

0.7 days 

1.4 days 

10 51 

r r = 4.7 days 

Td — 16 days 

2.8 days 

7.9 days 

1.1 days 
2.4 days 

0.5 days 

0.8 days 


the lifetime of the secondary star following mass accretion from 
the primary star, £2 ■ Given that the progenitor mass of the explod¬ 
ing ultra-stripped star is expected within the interval 8 — 25 Mg, 
we have that £2 — 7 — 40 Myr. (This is the original mass before 
this secondary star loses its hydrogen envelope in a CE, which cre¬ 
ates the naked helium star orbiting the NS, and which subsequently 
undergoes Case BB RLO.) 

For the prod uctio n of do uble neutron star (DNS) systems, it 
was argued by IVoss & Taurisl j2003h that the progenitor binaries 
must have been fairly wide following the first SN explosion to avoid 
a merger during the subsequent CE and spiral-in phase when the 
secondary star evolved. The reason is that only massive stars which 
are evolved (i.e. requiring a wide binary) have small enough en¬ 
velope binding energies to allow for a successful ejection of th eir 
envelope leading to surviva l of the CE-phase dDewi & Taurisl200ol . 
l200ll:lPodsiadlowskil200ll) . As a consequence of this criterion, also 
the pre-SN systems (prior to the first explosion) must, in general, 
have been wide. Therefore, to ensure survival of these systems, 
the first SN must have imparted a small kick on the newborn NS, 
less tha n 200 k ms -1 and with an average of only ~ 50 kms -1 
dVoss & Taurisl[2003l ). The resulting systemic velocities are there¬ 
fore also small with an average of v Bys — 10 kms -1 . This picture 
seems to have observational support when considering the typical 
small migration distanc es (~ 0.1 kpc) of FIMX Bs from the spiral 
arms in the Milky Way dColeiro & Chatvl[2013l) . 

To conclude, we find that the typical expected value of the 
travel distance is d = v sys • £2 — 200 — 300 pc (the entire inter¬ 
val spanning between 100 — 800 pc). Therefore, we predict ultra- 
stripped SNe to occur relatively close to the star-forming regions in 
host galaxies. 


4.3.1 Comparison to SN 2005ek 

In Paper I, we argued that SN 2005ek is a good candidate for a n 
ultra-stripped SN. This SN was discovered bv lDrout et al.l d2013l) . 
It exploded in the star-forming galaxy UGC 2526, but with a l arge 
projec ted offset to areas with strong Ha emission. IDrout et al.l 
d2013l) argued that, although the upper limit of the Ha line flux 
(3.3 x 10 37 ergs -1 ) is an order of magnitude below the mean 
value for Hll regions associated with core-collapse SNe dCrowtheH 
120131) . such a low value at the explosion site of core-collapse SNe 
has been measured previously. Furthermore, they argue that, if the 
progenitor star was older than ~ 20 Myr, the lack of observed 
Ha emission could be caused by a shorter lifetime for giant Hll 
regions. Given that our DNS modelling predicts progenitor life¬ 
times up to 40 Myr for ultra-stripped exploding stars (especially 
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Initial He —star mass, M Hei (M Q ) 

Figure 11. The maximum mass-transfer rate, |Mg^ x |, from the helium 
star donor during Case BB (or Case BA/BC) RLO as a function of its initial 
mass, MH e ,i, for binaries with different initial orbital periods (coloured 
lines). The light grey-coloured lines represent Kelvin-Helmholz mass- 
transfer rates, Mkh = ^He/ r KH< where tkh is the Kelvin-Helmholz 
(thermal) timescale of the donor star at the onset of the RLO. It is notable 
that iMjgf*! ~ 10 -4 Mq yr -1 for almost all values of Mn e ,i and dif¬ 
ferent initial orbital periods, and also that these rates are fairly close to the 
values of Mkh- The Eddington accretion rate of the NS is shown as a 
dashed line, some 3^t orders of magnitude smaller than |Mg^ x |. See text. 



Figure 12. The mass-transfer rate, |Mne| from the four helium star bi¬ 
naries discussed in Sections OElEUaiidlll] leading to an Fe CCSN, 
EC SN, ONeMg WD and a CO WD, respectively. The mass-transfer rates 
all peak close to 10“ 4 Mq yr —1 . The narrow peaks in \Mue\ at the 
end of the computations are caused by numerical instabilities and/or shell 
flashes. For the helium star with an initial mass of Mn e ,i = 2.6 Mq and 
P orb j = 0.06 days, the system leads to a merger in the UCXB-phase. The 
helium star with an initial mass of Mn e i = 3-5 Mq has reduced its mass 
to 3.0 Mq at the onset of RLO as a result of stellar wind mass loss (Fig. [3) 
For the other (low-mass) helium stars, stellar wind mass loss was neglected. 


for EC SNe), we find that SN 2005ek remains a good candidate for 
an observed ultra-stripped SN. 


5 HELIUM STAR MASS-TRANSFER AND NS 
ACCRETION RATES 

The calculated mass-transfer rates from our models are plotted in 
Fig. no This plot shows how the peak value of the mass-transfer 
rate, | Mh* x | during stable RLO (Case BA, BB or BC) depends on 
P or b,i and Mue,i of the helium star. For the smallest value of P 0 rb,i 
(Case BA RLO). |A/^ e ax | increases by more than an order of mag¬ 
nitude when Mne,i is increased from 2.5 to 3.5 Mq. However, for 
the wider systems with initial P or b,i ^ 2.0 days (Case BC RLO), 
AfJ I " ax is roughly constant or slightly decreasing with A-fHe.i- In 
general, |AT^f e ax | is always close to 10~ 4 Mq yr -1 within a factor 
of a few. We have also plotted the Kelvin-Helmholtz mass-transfer 
rate, Mkh, defined as the ratio between the stellar mass and the 
Kelvin-Helmholtz (thermal) timescale of the star at the onset of the 
RLO, i.e. Mkh = Mhb/tkh = RL/(GMu e ), where R and L are 
the radius and the luminosity of the star at the onset of the RLO, 
and G is the gravitational constant. The range of values of | M/J 1 ” | 
is fairly consistent with that of A/kh although the latter range has 
a larger spread. Only for very evolved helium stars in wide binaries 
with P 0 rb,i = 20 days is the value of AFkh — 10~ 2 Mq yr -1 (not 
shown on plot) significantly larger than | Afg 1 e ax |. However, whereas 
Mkh increases with P or b,i, the calculated values of |M^ rl e ax show 
the opposite trend for larger values of Mne.i- One reason for this 
is that the value of |Mg 1 e ax | is influenced by the shrinking of the 
orbit when mass is transferred from the (at least initially) more 
massive donor star to the lighter accreting NS, thereby leading to 
sub-thermal timescale mass transfer. 

Fig. ED shows examples of more detailed mass-transfer rates 
for the four systems discussed in Section [3 In general, the com¬ 
puted mass-transfer rates are close to 10 -4 M© yr -1 . This value 



Figure 13. The mass-transfer rate, \Mne I as a function of time for 3.0 Mq 
helium star donors with different values of P orb i- The time t = 0 is defined 
at RLO onset for all systems. For the two narrow systems the calculations 
are terminated by a vigorous helium shell flash following off-centre oxygen 
ignition, similar to the case discussed in Paper I. It is clearly seen how the 
duration of the mass-transfer phase, A t x , increases with decreasing value of 
Porb i> resulting in more rapidly spinning recycled pulsars in DNS systems 
in close orbits. In all cases shown the outcome is an ultra-stripped Fe CCSN 
with M enVj f = 0.04 — 0.23 Mq. For wider systems with P or b,i = 5 — 
120 days, see Fig. ED 

is 3^1 orders of magnitude larger than the Eddington accretion rate 
for NSs and therefore implies that > 99.9 per cent of the trans¬ 
ferred material is lost from the system (presumably in the form 
of a jet, or a wind from the disk, with the specific orbital angu¬ 
lar momentum of the NS, as assumed in the isotropic re-emission 
model applied to our computations). Depending on the exact as- 
su med accretio n efficiency for such Case BB RLO calculations (see 
[Lazarus et al ■ 20li for discussions), the NS thus typically accretes 
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~ 10 3 Mq of material (see also Fig. ED- For the values quoted in 
Table [I] we assumed an accretion efficiency of 100 per cent of the 
Eddington accretion rate. However, it is possible that the amount 
of material accreted by the NS, AMns, is larger by at least a fac¬ 
tor of three, based on evidence_for a high accretion efficiency from 
PSR J1952+2630 ^Lazarus et all2014h . 

Super-Eddington accretion by a factor of a few is eas- 
ily obtained by a thick accretion disk with a central funnel 
dAbramowicz et alii 1 98(1) . or a reduction in the electron scattering 
cross-section ca used by a strong B-field dBasko & Sunvaevll 19761 : 
IPaczvnskl 19921) . 


6 DOUBLE NEUTRON STAR (DNS) SYSTEMS 

6.1 Correlation between P or b and P spin for mildly recycled 
pulsars in DNS systems? 

In Fig.[l3] one can see from our computations how the duration of 
the mass-transfer phase, At x , increases as a function of decreasing 
Porb.i for a given binary (see also Fig. |T6] for wider systems). In 
other words, the larger P or b prior to Case BB RLO, the shorter is 
the duration of the mass-transfer phase (due to the more evolved 
stage of the donor star at the onset of the RLO) and thus the less 
efficient is the recycling process and thus the longer is the expected 
spin period, P sp m, of the mildly recycled pulsar. Assuming this re¬ 
lation to survive the SN explosion, one might expect a correlation 
between P or b and P sp m of the recycled pulsars in DNS systems 
(see further discussion in Tauris et al., in preparation). A correla¬ 
tion between P or b and P sp i n has previo usly been hinted indirectly 
on similar quali tative arguments (e.g. IChaurasia & Bailed 120051 : 
iDewi et all 20051) . However, we caution that the value of P ap i n de¬ 
pends on the location of the magnetospheric bound ary of the ac¬ 
creting NS and therefor e it depends on its B-field dPringle & Reesl 
1 1 9721 : [Lamb et al.fl973l) . This dependence is difficult to correct for 
given the current poor understanding of B-field decay via accretion 
and — probably worse — convolved with the scatter in NS B-fields 
from their birth and the unknown magnetic inclination angle which 
affects the braking torque and thus the observed time derivative of 

p . 

- 1 spin- 


6.1.1 Spin periods of mildly recycled pulsars in DNS systems 

As demonstrated in lTauris et alj ( 120121') . the calculated amount of 
material accreted on to NSs during Case BB RLO in post-CE bina¬ 
ries matches well with the observed spin period values for mildly 
recycled pulsars. Using their equation (14), we find that the values 
of AMns calculated here (between 10 -4 Mq and 3 x 10 -3 Mq, 
if assuming an accretion efficiency of 100 per cent up to the Ed¬ 
dington limit, cf. discussion in Section [5} are adequate to produce 
a mildly recycled pulsar with a spin period roughly between 25- 
350 ms. If the accretion e fficiency can be^super-Eddington by a 
factor of 1.5, 2.0 or 3.0 (Laza rus et al. I B . the fastest spin pe¬ 
riod of a DNS pulsar becomes 19, 15 and 11 ms, respectively. For 
comparison we note that the DNS system with the fastest spinning 
rec ycled NS known to date is the double pulsar PSR J0737—3039 
jBurgav et al.ll2003h which has a spin period of P sp j n = 23 ms. 

The formation and properties of very wide-orbit DNS systems 
(with post-SN P or b > 20 days) are discussed in Section ROl 



Final core mass, M coref (M Q ) 

Figure 14. Cross sections of pre-SN stars versus their final core mass. The 
plot is based on evolved helium stars with initial masses of 2.8 — 6.0 Mq 
and an initial orbital period of 0.1 days (cf. Table [7}. The total mass is 
divided into three regions representing the amount of helium (orange), 
the remaining carbon-rich envelope (red), and the metal core (purple). 
The blue shaded region represents stars leaving ONeMg WD remnants. 
The yellow line indicates the gravitational masses of the NS remnants 
(~ 1.20 — 1.80 Mq) calculated for a SN mass cut at the metal core bound¬ 
ary, ap plying the gravitational binding energy following lLattimer & Yahij 
<19891) . The hollow yellow hatching indicates that lower NS masses are 
possible if the mass cut occurs deeper in the exploding stars. Helium stars 
evolved in X-ray binaries with shorter initial orbital periods will be further 
stripped, resulting in even smaller envelopes (Figs. l9land[T0l. 


6.2 Neutron star masses and kicks 

In Fig. |T4l we present the interior composition of our pre-SN mod¬ 
els as cross sections by mass, which were calculated for helium 
stars with initial masses between 2.5 and 6.0 Mq, in binaries with 
an initial orbital period of 0.1 days and an accreting 1.35 Mq NS. 
(Note, the helium stars with the smallest initial masses do not pro¬ 
duce a SN, cf. Table [T}. This figure illustrates the extreme degree 
of stripping for, in particular, the lowest mass helium stars. We es¬ 
timate that the NSs produced from our ultra-stripped SNe could 
potentially have masses in the range 1.10 — 1.80 Mq. The ex¬ 
act mass of the NS depends on the mass cut during the explosion 
and the NS equation of state, both of which are still uncertain. We 
notice that the second formed NS in a DN S sy stem with the lowest 
known mass is PSR J0453+1559 dMartinez et al.l201^1 which has a 
precisely measured mass of 1.17 Mq (less than ~ 1.18 Mq at the 
99.7 per cent confidence level). For recent disc ussio ns on th e uppe r 
limit of NS birth masses in b i naries , see e.g. Ta uris e t al. d201ll) : 
lOzel et aljFoTinlFrver etaP ( l2012h . 

As a result of the SN explosion, a significant momentum kick 
is often imparted on a newborn NS. Such kicks are associated 
with explosion asymmetries and may arise from non-radial hydro- 
dynamic instabilities in the collapsing stellar core (i.e. neutrino- 


Blondin & Mezzacanna 2006!: 

Scheck et alj 120061; Foglizzo et al. 

2007 

; Marek & Janka 2003; 

Janka 2012|. 2013): Bruenn et al. 


2014 ; iFogfizz^^niSSF These instabilities lead to large-scale 


anisotropies of the innermost SN ejecta, which interact gravitation¬ 
ally with the proto-NS and accelerate the nascent NS on a timesca le 
of several seconds le.g. |jank3l2012l : lWongwathanarat et al.ll2013l) . 

The momentum kick imparted to a newborn NS via an EC SN 
is expected to be small. This follows from detailed simulations 
which imply: i) explosion energies significantly smaller than those 
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Figure 15. Gravitational binding energy of the envelope as a function of CO 
core mass (upper panel) and ONeMg core mass (lower panel) of our pre- 
SN models. Open and filled symbols correspond to models leading to an 
EC SN and an Fe CCSN, respectively (cf. Section [3}. The different colours 
represent P 0 rb,i °f the original helium star. The error bars were calculated 
by varying the core/envelope boundary by 0.01 Mq. The black arrow in¬ 
dicates the change in location in the diagram for one particular model, from 
the stage of early core oxygen burning to a more evolved stage near the end 
of core oxygen burning. Two data points (shown with arrows at the top of 
the lower panel) are located at Pbind - 1-4 x 10 50 erg. 


inferr ed for standard Fe CCSNe jKitaura et al.ll2006l ; lDessart et al.l 
l2006t) : and ii) short timescales to revive the stalled supernova 
shock, compared to the timescales of the non-radial hydrodynamic 
instabilities which are requ i red to produc e strong anisotropies, e.g. 

Simulations of EC SNe 


instabilities which are required to produce 
IPodsiadlowski et all ]2004l):ljankal d2012h. 


J 2 * 00 ftl ) predict explosion t 


by Kitaura et alj 1 2006h : Dessart et al 
ergies of about (or even less than) 10 50 erg, and thus possibly kick 
velocities less than 50 kms -1 , i.e. significantly smaller than the 
average kick velocities of the order 400 — 500 k ms -1 which are 
evidently imparted on y oung pulsars in general jLvne & Lorimetl 
11994 iHobbs et alj|2005l) . 


One may ask what would be the magnitude of a kick imparted 
to a NS born in an ultra-stripped SN, which can be either an EC SN 
or an Fe CCSN. We see two factors which imply that in these ob¬ 
jects, the NS kicks may by small. First, from our modelling of the 
progenitor stars of ultra-stripped SNe, we have demonstrated (Sec¬ 
tion 0 that the amount of ejecta is extremely small compared to 
standard SN explosions. This may likely lead to a weaker gravi¬ 
tational tug on the proto-NS (caused by asphericity of the ejecta, 


e.g. lJankall20l3) and thus a small kick. Second, the binding ener¬ 
gies of the envelopes of our models, even including the CO-layer, 
are often only a few 10 49 erg (Fig. ED, such that even a weak 
outgoing shock can quickly lead to their ejection, potentially be¬ 
fore large anisotropies can build up. While these arguments hold 
best for the lowest mass NSs formed from our models, it is pos¬ 
sible in principle that one can form even rather massive NSs from 
ultra-stripped progenitors with low kicks, provided that the whole 
(most) of the ONeMg core collapses. However, these more massive 
ONeMg cores are also expected to produce larger iron cores and 
may therefore lead to more ’normal’ Fe CCSNe with larger explo¬ 
sion energies and larger kicks. We are unable to resolve this issue 
as this requires full core-collapse calculations. In addition, it will 
be left to a future study to investigate how the final pre-SN cores in 
ultra-stripped progenitors differ from those of single stars of com¬ 
parable initial mass to confirm whether low-kick NSs preferentially 
occur in close binaries, as seems to be indi cated by observations 
( iPfahl et alj 2002 i :FP odsiadlow ski et al .112004 ). 

While it is difficult to estimate kick values, these argu¬ 
ments taken together allow us to speculate that ultra-stripped SNe 
(EC SNe and, at least, Fe CCSNe with relatively small metal cores) 
generally lead to the formation of NSs with small kicks. This would 
have important consequences for the formation of DNS systems 
and we discuss these implications in a future study. 

We recall that the metal core masses (M corei f ) listed in TableQ] 
correspond to the mass inside the CO core boundary. Whereas this 
boundary is fairly stable towards the end of the nuclear evolution of 
our modelling, the boundary of the ONeMg core is quite sensitive 
to the number of carbon burning shells which develop in the final 
phases of our calculations. This explains the somewhat larger scat¬ 
tering in envelope boundary energies in the lower panel compared 
to the upper panel of Fig. ED 

6.3 Formation and properties of very wide-orbit DNSs 

In the previous sections, we have discussed the evolution of tight 
orbit (P 0 rb,i ^ 2 days) helium star-NS systems, leading to the for¬ 
mation of fairly close orbit DNS systems via an ultra-stripped SN. 
In addition, we have investigated a number of wide-orbit helium 
star-NS systems with 5 ^ -Porb.i ^ 120 days, cf. Table]]] 

A common feature of these wide-orbit systems is that by the 
time the helium star fills its Roche lobe (at an orbital period larger 
than P 0 rb,i as a consequence of stellar wind mass loss), the helium 
star is quite evolved and is already undergoing carbon shell burning. 
As a consequence of this, the remaining lifetime of the helium star 
is less than a few thousand years at the onset of RLO. Therefore, 
by the time the helium star explodes in a SN (during RLO) it will 
still have a rather thick helium envelope of = 0.2 — 0.8 Mq 

(the wider the orbit, and the more massive the initial helium star, 
the larger envelope will remain at the onset of the core collapse). 
Such a relatively massive envelope ejected at the time of the SN 
will not only affect the observational properties of the light curve. 
It also results in larger post-SN orbital periods and eccentricities 
of the DNS systems formed. Typical values for the eccentricities in 
the case of symmetric SNe would be 0.2 — 0.5. However, larger 
(and smaller) values are possible if kicks are added (which indeed 
might be the case as a result of the larger ejecta mass of these stars, 
cf. discussion in Section 1672b . We find that orbital periods of DNS 
systems with a recycled pulsar are expected out to about 200 days 
(again depending on the imparted kicks). 

Given that the SN will terminate the recycling process of the 
accreting NS shortly after RLO is initiated, we expect the resulting 
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Initial orbital period, P orb; (days) 

Figure 16. Remaining helium envelope mass, (red curve) at the 

time of the SN as a function of initial orbital period, P or b i of helium 
star-NS systems with Mn e ,i = 3.0 Mq. Also plotted is the duration of 
the RLO prior to the explosion (blue curve), which determines the proper¬ 
ties of the (marginally) recycled pulsar companion (see text). In all cases 
shown, the metal core reaches a final mass of ~ 1.58 Mq at the mo¬ 
ment of the Fe CCSN. Ultra-stripped SNe require P orb i 2.0 days for 
Mjje i = 3.0 Mq. The grey-shaded region to the right indicates the loca¬ 
tion of systems which are too wide for post-HMXB/CE Case BB RLO to 
occur. 


wide-orbit DNS systems to contain extremely mildly recycled pul¬ 
sars, i.e. marginally recycled pulsars, with spin periods in excess 
of 100 ms and relatively large B-fields. These pulsars should be 
observed in the grey-zone in the PP-diagram where it is difficult 
to distinguish such marginally recycled pulsars from the lower El- 
field population of isolated, old non-recycled pul sars. The recently 
discovered DNS system PSR J1930—1852 iSwiggum et aljfiola) 
is an bona fide example of such a system. It has a spin period of 
185 ms, a spin period derivative of P = 1.8 x 1CD 17 ss -1 , an 
orbital period of 45 days and an eccentricity of 0.40. 

In Fig. [16] we have plotted the amount of helium left in the 
envelope at the time of the SN, as well as the duration of the 
pulsar recycling phase via Case BB (strictly speaking Case BC) 
RLO, as a function of the initial orbital period of the helium star- 
NS system. In this plot, calculations are shown for helium stars 
with Mue,i = 3.0 Mq - see Table |T| for other systems with 
P orb ; = 20 days. It is clearly seen how the remaining helium en¬ 
velope mass increases with P or b,i, and how the duration of the RLO 
(Af x ) decreases with P or b,i- Given that the growth of NS mass is 
expected to be limited by (possibly a few times) the Eddington- 
accretion limit of MEdd.NS ~ 4 x 10 -8 Mgyr -1 (for pure he¬ 
lium), we expect wide-orbit DNS systems to be (marginally) recy¬ 
cled with AMns — 1 — 10 x 10 -4 Mq, typ ically resulting in spin 
periods in excess of 100 ms [cf. eq.(14) in lTauris et al.l J20 1 21) for a 
rough estimate of the expected spin period after recycling]. 

The system with P or b,i = 120 days will not initiate RLO 
since the radius of the helium star, reaching a maximum of 74 Rq, 
is unable to catch up with its Roche-lobe radius (which is increasing 
since the orbit is expanding due to stellar wind mass loss). Hence, 
its NS companion will avoid recycling. 

Finally, we note that in all cases, even for the widest sys¬ 
tems which initiate RLO, the expanding helium envelope is stable 
against convection (except in a very tiny layer near the surface). 
This means that: i) the mass-transfer rates in these systems are not 
much different from those shown in Fig. |TT] for helium star donors 


in closer orbits, and ii) the RLO is always dynamically stable such 
that these systems avoid evolving into yet another CE. 


6.4 Constraints on NS accretion prior to Case BB RLO 

Following the formation of the first born NS in a massive binary, 
there are four phases of evolution which could potentially lead to 
accretion onto this NS: i) accretion from the stellar wind of the mas¬ 
sive donor star prior to the formation of a CE (i.e. standard HMXB 
accretion), ii) accretion during the CE-phase, iii) accretion from the 
stellar wind of the exposed naked helium core (Wolf-Rayet star) of 
the former massive companion star, and iv) Case BB RLO from the 
expanding helium star as investigated in this paper. 

Depending on the lifetime (mass) of the companion star, an 
upper limit to the amount of material accreted from a wind in the 
first phase is of the order 10 -2 Mq. The wind accretion from the 
naked helium star prior to Case BB RLO (third phase) may reach 
a few 10~ 4 Mq. However, in both cases it is unclear to which 
extent the wind-accreting NS is able to efficiently spin up. This de¬ 
pends on, for example, whether or not an accretion dis k is f ormed, 
and w hether or not propeller effects are important (e.g. lFrank et al.l 
120021 ) . 

It was recently demonstrated by iMacLeod & Ramirez-Ruizl 
d2015allbl) that a NS embedded in a CE only accretes a modest 
amount of material during its spiral-in as a result of density gra¬ 
dients across the NS’s accretion radius (which strongly limits ac¬ 
cretion by imposing a net angular momentum to the flow around 
the NS). They computed a firm upper limit of 0.1 Mq. Also in 
this phase, given the extreme conditions with a very high ram pres¬ 
sure and a squeezed magnetosphere, the efficiency of transporting 
angular momentum to the NS is likely to be very limited. 

We therefore conclude that the properties of the first 
born NS in an observed DNS system might entirely be ex¬ 
plained by the Case BB RLO investigated here. Further evi¬ 
dence for this comes from observation s of the wide-orbit DNS 
systems: PSR J17 53-2240 dKeith et alJliOQ^). PSR J1811—173 6 
dLvne et alj|200ch and PSR J1930-1852 dSwiggum et~ahll2015l) . 
which have rather slow spin periods of 95, 144 and 185 ms, respec¬ 
tively, in agreement with our calculated estimates from Case BB 
RLO and indicating rather inefficient angular momentum transport 
during the other phases briefly discussed above. 


6.5 Merger time of DNS systems resulting from our models 


As previously mentioned, the final orbital period prior to an ultra- 
stripped SN, P or b,f, is typically between 1 hr and 2 days (cf. Ta¬ 
ble |T]l. In general, the dyna mical effect s of S Ne may chang e the 
orbital period significantly dHillsI 19831: Brandt & Podsiadlowskil 


1 19951 : iTauris & Takendn~998l : IWex et alJl200(lh . However, the core 

structure of a star depends strongly on whether it evolves in a c lose 
binary or in a wide binary/isolation dPodsiadlowski et al.l2004i ). As 
suggested in Section [6Al kicks imparted on NSs formed from ultra- 
stripped stars (i.e. during the second explosion forming DNS sys¬ 
tems) are possibly small. If correct, this would imply that, on av¬ 
erage, the post-SN orbital periods will not be much larger than the 

pre-SN orbital periods. _ 

As a consequence of gravitational wave radiation J Peter. si 
Il964h . the merger timescale for such newly formed DNS systems 
with orbital periods between 1 hr and 2 days is between a few Myr 
and oo; shortest for close binaries with high post-SN eccentrici¬ 
ties (which may still be significant even from small kicks). The 
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Merger time, log r gwr (Myr) 

Figure 17. Histogram of merger times caused by orbital decay via gravita¬ 
tional wave radiation of the final double degenerate systems resulting from 
the binaries modelled in this paper with P Q rb i = 0.1, 0.5 and 2.0 days (cf. 
TableQ), corresponding to the three peaks, respectively. The blue colour 
represents DNS systems; the red colour represents NSWD systems (i.e. 
systems where the helium donor star left behind a massive WD orbiting 
the NS). See text for details. 



Figure 18. Grid showing the outcome (the final fate of the helium star 
donors) of our modelling as a function of the initial values of orbital period 
and helium star mass - see also Table Q] The dotted lines (equations 13- 
Go) separating the regions of different outcomes are useful for population 
synthesis. Our modelled systems with P or b ; of 5 — 120 days are not 
shown here, but fall in the region between P or b,i = 2.0 days and the 
single helium stars. 


maximum post-SN orbital period of a DNS system that will merge 
within a Hubble time (assuming NS masses of 1.4 M 0 ) is between 
0.7 and 1.7 days for post-SN eccentricities of 0.1 — 0.7. 

In Fig. [TT] we have plotted the merger times, r gwr , caused 
by orbital decay via gravitational wave radiation of the double 
degenerate systems resulting from our models presented in Ta¬ 
ble Q] In case the donor star terminated as a SN, creating a DNS 
system, we assumed, for simplicity, a symmetric explosion (i.e. 
without a kick). The mass of the newborn NS was calculated us¬ 
ing: Mns = A'fcore.f — where the released binding en¬ 

ergy of the NS correspon ds to a mass decrease of AT^g 1 ' 1 = 
0.084 M 0 • (Mns/ATq ) 2 iLattimer & Yahill989ll . i.e. following 
the yellow line of NS masses in Fig. El Subsequently, we cal¬ 
culated the amount of mass eject ed a nd the result in g post-SN or¬ 
b ital period and eccentricity (e.g. iBhattacharva & van den Heuvell 
1991), before calculating T gwr via numerical integration IPetersI 
1964). In case the donor star terminated as a WD, we assumed 
that Mwd = AT core f and that the orbit remained circular with 
a resulting orbital period equal to P or b,f. (The slight widening of 
the orbit from losing the thin envelope in a proto-WD wind is as¬ 
sumed to be compensated with orbital decay due to gravitational 
wave radiation during this epoch.) The systems in the histogram are 
seen to cluster in three groups corresponding to r gwr < 250 Myr, 
r gwr = 1 — 10 Gyr and r gwr > 85 Gyr, reflecting our sys¬ 
tems chosen initially with helium-star ZAMS orbital periods of 
Porb,i 4 0.1 days (left), 0.5 days (middle), and 2.0 days (right), cf. 
Table E] We emphasize that the distribution in Fig. Ell only repre¬ 
sents examples within the expected range of r gwr and that the true 
dist ri butio n o f r gwr requires population synthesis modelling (e.g. 
IVoss & Taurisll20()3l) . It is possible that some of the tight systems 
which end up in a CE may survive the spiral-in phas e and produce 
very tight systems with even smaller values of r gwr l iDewi & Polsl 
120031) . see discussion in Section [7~2l 

In a future publication we analyse in more detail the observed 
distributions of orbital periods, spin periods, eccentricities and sys¬ 
temic velocities of DNS systems and confront them with simula¬ 
tions and expectations from ultra-stripped SNe. 


7 DISCUSSIONS 

7.1 Mapping initial parameter space to final outcome 

Fig. ED shows the final fate of the 2.5 — 3.5 M 0 helium stars 
which we have modelled in binaries. Also shown to the right is 
the outcome of the same helium stars evolved in isolation without 
a companion star. As discussed in Section 13.11 the boundary be¬ 
tween EC SNe and Fe CCSNe is perhaps somewhat less certain 
compared to the other boundaries shown. Furthermore, along the 
sequence calculated for P OI b,i = 0.06 days we could not deter¬ 
mine (due to numerical instabilities) the exact mass for which the 
mass transfer becomes dynamically unstable, leading to a CE evo¬ 
lution. Certainly, the helium stars with M-a e y > 3 M 0 end up in 
a CE for these short initial orbital periods. Systems with P or b,i less 
than 0.06 days are excluded because the helium ZAMS stars cannot 
fit inside such small Roche lobes. One can use this figure for deter¬ 
mining the fate of helium stars in close binaries when performing 
population synthesis. 

7.1.1 Fitting formulae with applications to population synthesis 

Based on our results presented in Figs.[2][9]En] an dEH] we have fitted 
a number of relations which are, for example, useful for population 
synthesis modelling. 

For the relation in Fig. [2] between final pre-SN metal core 
mass, M COIB} f and initial values of helium star mass, A/He,i (prior 
to Case BB RLO and wind mass loss), and orbital period, P or b,i we 
find: 

“ (srb- +a49 ) (j£ " 01 “)' (3) 

where P or b,i is in units of days and masses are in Mq. For the 
investigated mass interval of 2.5 < Afn e ,i/AT 0 ^ 3.5 this fit is 
accurate to within about 1 per cent (except for the very low-mass 
end where the discrepancy is 3 per cent for the wider systems). 

The remaining amount of helium in the pre-SN envelope, 
can be fitted as a function of final pre-SN metal core mass 
and initial orbital period as follows: 
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wenv ^ f 0 - 18 P™ ■ [HMU f) - 1.05] , 

He ’ f \ M core ,f • [ln(P o 7 b °. 2 ) + l] + ln(P 0 rb,i) - 1-5, 

where the upper expression is valid (generally within an error of 
±0.01 Mq) for 0.06 < Porb.i ^ 2.0 days, and the lower expres¬ 
sion is an approximation for P or b,i > 2.0 days and single helium 
stars (P or b,i —> oo). By combining equations 0 and 0 one can 
calculate directly as a function of initial orbital period and 

initial helium star mass. 

The orbital evolution can easily be calculated analytically 
by separating the wind mass loss prior to Case BB RLO from 
the subsequent Roche-lobe mass transfer (during which the effect 
of wind mass loss is small compared to a mass-transfer rate of 
~ 10 -4 Mgyr -1 , cf. Section 0. For mass loss due to a jiirect 
fast stellar wind from the h elium star (i.e. Jeans’s mode, iHuand 
1 1963l : lvan den Heuvelll 1 9941) . the orbital period change is given by: 


P_ = (l + q o \ 2 
Po \ i + q ) 


(5) 


where Po and P represent the initial helium ZAMS orbital period 
(Porb.i) and the orbital period at the onset of the RLO, respectively. 
The mass ratio between the stellar components, q = Mhb/Mns at 
the same two epochs is denoted by qo and q. 

By integrating the orbital angular momentum balance equa¬ 
tion and applying th e so -called isotropic re-emission model 
jBhattacharva & van den Heuvelll 199lh one can find the change in 
orbital separation during non-conservative RLO. lThurisI dl9960 de¬ 
rived: 

3g~5 r. 

a _ ( qo(l ~ 13) + l\^ ir ( qo + 1\ ( qoY , 

a 0 U(i-/3) + i; U + i/ \q) ’ 


where ao and a refer to the orbital separation prior to RLO and 
after RLO, respectively, and where qo and q now represent the 
mass ratios at these two epochs. The parameter /3 represents the 
(constant) fraction of transferred material which is lost from the 
vicinity of the NS. Given that the mass-transfer rate of Case BB 
RLO is much higher than the Eddington limit for an accreting NS 
(IM 2 I AfEdd,NS> an d thus f3 —> 1) one can approximate the 
abo ve equation and c ombine it with Kepler’s third law to yield (see 
also lKing et al .1120(771) : 


lim 

/3— >1 





(7) 


As a numerical example, we can compare this expression with our 
BEC computation of Case BB RLO for a 3.2 Mq helium star orbit¬ 
ing a 1.35 Mq NS with an initial orbital period of 2.0 days, and for 
which stellar winds were not included. As seen in Table [T] the final 
mass of the donor star prior to the explosion is 1.98 Mq . Hence, we 
get qo = 2.37037 and q = 1.46656 for this system, and inserting 
Po = 2.0 days in the above equation we find that the final orbital 
period at the end of Case BB RLO is P = 1.0476 days (exactly 
as calculated with our binary stellar evolution code, cf. 1.05 days 
in Table |T}. Of course, gravitational wave radiation included in the 
code also carries away orbital angular momentum. However, this 
effect is negligible during the short lasting (typically a few 10 4 yr) 
Case BB RLO, cf. Fig. [13] 

For the final fate of a helium star in a close orbit with a NS 
(Fig. 00 dotted lines) we find that the outcome is either: 


Fe CCSN : 

EC SN : 
ONeMg WD : 


Mne,i > MFeCCSN, 

MeCSN ^ MHe.i < MFeCCSN, 
MHe.i < MeCSN, 


( 8 ) 


where the threshold initial mass for producing an Fe CCSN is: 

MFeCCSN = 2.74 P ~™° 209/Polh ' 1 , (9) 

and the threshold initial mass for producing an EC SN is: 

Mecsn = -r4-+ 2 - 58 - (10) 

61 Porb.i 

where P or b,i is in days and the threshold masses are in Mq. The 
above formulae are valid for P or b,i > 0.07 days and MHe,i ^ 
2.5 Mq. For P or b,i — 0.06 — 0.07 days, the outcome is a CO WD 
if Mn e ,i 2.8 Mq, and most likely a CE for more massive he¬ 
lium stars. Furthermore, one should keep in mind that all calcula¬ 
tions were performed for a 1.35 Mq accreting NS. Using other NS 
masses may slightly change the boundaries stated above. 


7.2 Comparison to previous work 


Whereas ultra-stripped SNe as such have not been discussed much 
in the literature up to now (besides from Paper I), we can compare 

our Case BB RLO modelling to earlier work._ _ 

The sy s temat ic in vestigations by _ Dewi et aTI d2002h . 

iDewi & Polsl d2003h and llvanova et al.1 1 12003 ). also considered 


Case BB RLO from helium stars with masses of 1.5 — 6.7 Mq in 
binaries with a NS. However, whereas these previous computations 
usually terminated their calculations after carbon burning, we were 
often able to continue the binary stellar evolution until oxygen 
burning. Hence, we have been able to strip the envelopes of our 
donor stars further prior to the SN (i.e. ultra-stripped SNe), in 

some cases all the way down t o < 0.01 Mq ._ 

The main conclusion of IDewi & Polsl d2003l ) is that helium 
stars with Ma e ,i = 2.8 — 3.3 Mq, and those with MHe,i = 
3.3 — 3.8 Mq in narrow systems of P or b,i < 0.25 days, lead 
to unstable RLO and formation of a CE. They argued that the 
outcome of these systems is ejection of the helium envelope, and 
thus survival of such a (second) CE-phase, and eventually forma¬ 
tion of tight orbit DNS systems with P or b ~ 0.01 days. For the 
more massive helium stars with MHe.i > 3.8 Mq, and those 
with MHe.i = 3.3 — 3.8 Mq in somewhat wider systems of 
Porb.i > 0.25 days, they found stable RLO leading to final 
DNS systems with P or b ~ 0.1 — 1 days. The outcomes of our 
computations are quite different. From Fig. |T 8 ] and Table Q] it is 
clearly seen that all our low-mass (2.5 — 6.0 Mq) helium stars 
are dynamically stable throughout the RLO. The only exceptions 
are for MHe.i ^ 3.0 (2.8) Mq in very close initial orbits of 
fo r b,i = 0.06 days where mass transfer is initiated early on the 
main sequence (Case BA). In some cases our models experience 
numerical instabilities near the end of the RLO, for example, trig¬ 
gered by a vigorous helium shell flash after ignition of oxygen (see 
Paper I). However, as demonstrated here (cf. Section [7~6l and the 
Appendix), even in case such events lead to a runaway mass trans¬ 
fer and the possibility of forming a (second) CE, the remaining 
lifetime of the donor stars until the SN explosion in these systems 
is much shorter than the timescale of the spiral-in process in a di¬ 
lute envelope with a mass of typically < 0.01 Mq. In addition, 
the binding energy of such an envelope is very small, leading to 
rapid ejection without significant spiral-in. Therefore, we do not 
(directly) produce such extre mely short P OI b DNS sys tems as pre¬ 
dicted by IDewi et all J2002h and IDewi & Polsl d2003l) . Of course, 
all our DNS systems formed with P or b < 0.5 days will also be¬ 
come extremely tight systems within a few Gyr due to gravitational 
wave radiation (Fig. Q7>, so this issue should not affect much the 
predictions for the detection rate of LIGO/VIRGO which are based 
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on population synthesis (e.g. lBelczvnski et al.ll2002l : IVoss & Taurisl 
1200 3l) . However, this difference in characteristics for the newly 
formed DNS systems is quite important for the numbers of observ¬ 
able DNS systems as radio pulsars with P or b < 1 day. 


The results of llvanova et all J2003I) ar e somewhat c l oser to 
ours. A remaining key difference is that llvanova et all 1200 3lt 
find that their calculated mass-transfer rates for systems with 
initial orbital periods, Porb.i < 0.4 days sometime exceed a 
critical mass-transfer rate, M c r it related to th e location of the 
so-called trapping radius (e.g. Begelmanl 1 19791 : IChevalierl 1 19931 : 
iMacLeod & Ramirez-Ruizl l2015ah . In systems with very large 
super-Eddington mass-transfer rates, matter presumably piles up 
around the NS and forms a growing, bloated cloud engulfing a 
large fraction of the accretion disk. A system will only avoid a 
CE if it manages to evaporate the bulk of the transferred matter 
via the liberated accretion energy at a distance from the NS larger 
than the trapping radius. Otherwise, the incoming material has too 
much negative binding energy to be ejected. At the same time, this 
trapping radius must be located inside the Roch e lobe of the NS 
in order to avoid a CE ( [King & Begelmanll 1 9991) . The exact loca¬ 
tion of the trapping radius, and thus the value of M cr it, is diffi¬ 
cult to calculate beca use i t also d epends on th e coo ling proce sses 
of the infalling gas dNaravan & Yil 1 19951 : iBlandford & Begelmanl 
Il999ll . Nevertheless, in all of our models presented in Table[l]with 


Porb.i > 0.06 days we find that |Mh 7 x | ^ Merit, where Merit 
is calculated from equation 15 in llvanova et al.1 (20031) . and varies 
from about 3.7 x 10 -3 M 0 yr -1 (for P or b,i = 2.0 days) to 
4.3 x 10“ 4 Mq yr -1 (for P or b,i = 0.08 days). Only for P or b,i = 
0.06 days we find a couple of systems where | M ^ x | > M cr it. 
However, these systems are anyway found to result in a runaway 
mass transfer and thus formation of a CE. Therefore, our calculated 
mass-transfer rates from the hel ium stars dono rs must be slightly 
smaller than those calculated by llvanova et alj d2003l ). The reason 
for this, besides from using different stellar evolution codes, is per¬ 
haps t hat our mass- transfer r ates are calculated us ing the prescrip¬ 
tion bv lRitted dl988h whereasllva nova et ail (2003 ) adopted the pre- 
scription by iTout & Eggletoid j 198 8l) . To summarize, we caution 
that all numerical calculations of Case BB RLO could potentially 
be affected by uncertain accretion processes at high mass-transfer 
rates exceeding the Eddington limit by ~ 4 orders of magnitude. 


De wi et alj J2002h evolved a number of helium star donors to 
become CO WDs or ONeMg WDs. For Case BA RLO they found 
that all helium stars with Mu e ,i = 1.5 — 2.9 M 0 form CO WDs. 
This is in good agreement with our calculations. We find that he¬ 
lium stars with a mass up to (at least, possibly slightly higher) 
2.7 Mq will produce CO WDs with amass of 0.78 —0.87 M 0 .For 
Case BB RLO they find that 2.1 — 2.5 M 0 helium stars produce 
ONeMg WDs. This is also in fine accordance with our calculations 
where we find that ONeMg WDs (of masses 1.21 — 1.35 M 0 ) 
are produced from helium stars with an initial mass of up to 
2.6 - 2.7 Mq for Case BB RLO (cf. Fig.[l8]». 


To determine the outcome of our models at the borderline be¬ 
tween EC SNe and Fe CCSNe we can compare our comp utations 
to the work of lUmeda et al.l 1 201 2l i and I Jones et al.l J20131) . This is 
necessary at this stage since we could not calculate through silicon 
burning until the onset of the core collapse. For example, we no¬ 
tice that the computed structure of our 2.9 Mq helium star with 
forb.i = 0.1 days (discussed in detail in Paper I) resembles very 
well the core of 9.5 — 11 M 0 single stars computed by these au¬ 
thors and which were found to produce Fe CC SNe. 


7.3 Varying the mixing length and the metallicity 

In our default models, we assumed a mixing-length parameter of 
a = l/H p = 2.0 and there fore we tested a number of models 
using a = 1.5 (Langeilfl 99 ll ). For close-orbit systems (P or b,i = 
0.1 days), we did not find any differences, whereas for very wide 
systems (P or b,i = 20 days) the stripping of the helium envelope 
was more efficient when applying a = 1.5, resulting in reduced 
envelope masses of about 25 per cent. However, in these cases the 
final metal core masses were only slightly smaller by < 2 per cent, 
and thus the value of the mixing-length parameter is not very cru¬ 
cial for the results presented in this study. 

We computed models with a helium star metallicity of Z = 
0.02. Given the possibility for long delay times between formation 
and merger events, it is possible that a fair fraction of the merg¬ 
ing DNS systems that LIGO/VIRGO will detect were formed in an 
evolutionary phase of the local Universe where the metallicity was 
lower. A decrease in metallicity results in smaller (and bluer) stars 
with a weaker stellar wind mass loss. Hence, compared to more 
metal-rich environments, the outcome is smaller stars in systems 
which widen less due to stellar winds. These two competing effects 
can, in general, be mimicked for Case BB RLO by a more metal- 
rich star in a different initial orbit. 


7.4 The ZAMS progenitor masses of our naked helium stars 


The ZAMS mass interval of the original hydrogen-rich progenitor 
stars, of the naked helium stars studied here, is difficult to determine 
precisely since it depends on the initial ZAMS binary separation 
and the modelling of convective core-overshooting, rotation, chem¬ 
ical mixing, metallicity, as well as the physics of the ejection pro¬ 
cess of the CE. From a de nsity profile analysis of the stellar models 
evolved bv lBrott et al.l d201 If) , we find that a 10 M 0 ZAMS star 
(Z = 0.02, non-rotating) produces a helium core of 2.4 — 3.2 M 0 
at the tip of the red giant branch, for a convective core-overshooting 
parameter in the interval $ov = 0.1 — 0.335. Similarly, a 20 M 0 
ZAMS star produces a helium core in the interval 7.7 — 8.9 Mq. 
These helium cores will be exposed via the CE and spiral-in phase 
as soon as the progenitor star fills its Roche lobe in the HMXB 
system. The thickness of the remaining outer hydrogen-rich layer 
at this point is somewhat uncertain. The location of the hydrogen 
shell source in these progenitor stars is initially located at smaller 
mass coordinates and thus the final (pure) helium star masses are 
slightly smaller. Furthermore, these (now) naked helium cores, or 
Wolf-Rayet star s, will subsequently reduc e their masses further by 
wind mass loss l lWellstein & Langeill 19991) . 

The initial masses of our binary stars which produce 


(Podsiadlowski et al. 2004: 

Umeda et al. 20121: Jones et al. 

20131, 

20 ll; Takahashi et alf 2013 

; Dohertv et al. 20151 and references 


therein). It has previously been demonstrated that the initial 
mass rang e for produ cing EC SNe is smaller in single stars 
dPoelaren ds et al.l 2008l) compared t o binary stars in non-degenerate 


systems 1 Podsiadlowski et ajj2004l) . 


7.5 Stripping by a black hole or a WD companion star 

Replacing the accreting NS with an accreting black hole (BH) 
should also lead to ultra-stripped SNe in close binaries. As a conse¬ 
quence of their larger masses, we expect accreting BHs would give 
rise to slightly smaller mass-transfer rates (more stable Case BB 


© 0000 RAS, MNRAS 000, 000-000 
















































































18 Tauris, Langer & Podsiadlowski 


RLO) and also lead to wider pre-SN orbits. However, given the ini¬ 
tial mass function, we expect less ultra-stripped SNe to occur in 
binaries with BHs. This is also reflected in the hitherto missing de¬ 
tection of a radio pulsar-BH binary (the descendent system of such 
a helium star-BH binary). 

Whether or not a massive WD is able to strip a helium star ef¬ 
ficiently prior to a SN is more questionable. We k now of two WD - 
NS systems, PSR B2303+46 (P or b = 12 davs. lStokes et all 19851 : 
Ivan Kerkwiik & Kulkarnil 1 19991) and PSR J1141—6 545 (P orb = 
0.20 days, Kaspi et alj2000l : Antoniad is et alj201 1 where the ob ¬ 
serve d radio pulsar was created after the WD l lTauris & Sennelsl 
|2000l) . Hence, it is quite likely that these systems evolved through 
a post-CE Case BB RLO phase where a massive WD was indeed 
accreting from an evolved helium star prior to its explosion. It is 
also clear in these two cases that the systems remained dynamically 
stable and that the WDs were able to survive the presumably highly 
super-Eddington mass-transfer rate (~ 10 -4 Mg yr~ 4 , Section]?} 
from their helium star companion without merging. How the ex¬ 
cess material was ejected from the system in these cases is not well 
understood. However, it is possible that significant stripping of the 
helium star prior to the SN was at work too in both of these systems. 


7.6 Effects of a dilute envelope on the pre-SN evolution 

As also discussed in detail by iDewi & Polsl j2003l) and 
llvanova et alj l l2003l) . runaway mass transfer in helium star-NS bi¬ 
naries may lead to the onset of a (second) CE evolution. The final 
outcome is determined by the competition between the timescale 
for spiral-in of the NS and the remaining lifetime of the pre-SN 
donor star before gravitational collapse. In general, the duration 
of the CE and _spiral-in phase is found to be < 10 3 yr (e.g. 
IPodsiadlowskilPOOll) . which is long compared to the estimated fi¬ 
nal lifetime of our models (in some cases just a few days). In the 
Appendix we estimate the timescale of NS spiral-in caused by a 
helium shell flash following ignition of oxygen burning. We find 
that the spiral-in timescale is often a factor of (at least) 100 longer 
than the remaining lifetime of the evolved helium star until the SN 
explosion. Hence, it is safe to ignore the expanding envelope at the 
late stages of evolution when determining the final outcome of the 
models. 


8 SUMMARY 

We have performed a systematic study of the evolution of close bi¬ 
naries containing a helium star and a NS. These systems are the 
descendants of HMXBs which undergo CE and spiral-in evolution 
and expose the core of an OB-star as a naked helium star. The sub¬ 
sequent evolution of these helium stars causes an additional mass- 
transfer towards the NS, often during helium shell burning (i.e. 
so-called Case BB RLO) when they expand to become giant stars 
again. In particular, we have investigated the expected nature of the 
resultant SN explosion (EC SN vs Fe CCSN) as well as the amount 
of stripping of the helium star donors prior to the explosion in order 
to predict the observational properties of these events as SN Ic or 
SNIb. 

We studied helium stars with masses between 2.5 — 10.0 Mq 
with an emphasis on the lighter ones with masses of 2.5 — 3.5 Mq 
which expand more during their evolution. In all cases we assumed 
an accreting NS with an initial mass of 1.35 Mq and varied the 
initial orbital period between 0.06 — 120 days, leading to either 
Case BA, Case BB or Case BC RLO. For comparison, we also 


evolved a sequence of isolated helium stars. The final fates of our 
helium stars span from CO WDs and ONeMg WDs to NSs pro¬ 
duced via EC SNe and Fe CCSNe. In the latter case, we were able 
to compute the evolution all the way to the onset of silicon burn¬ 
ing, a few days prior to the gravitational collapse. The results of 
our study are presented in a table (Table Q} with several parameters 
and a chart (Fig. 118b with the mapping of initial parameter space 
to final outcome. In addition, fitting formulae are derived with easy 
applications to, for example, future population synthesis studies. 

We define ultra-stripped SNe as exploding stars which contain 
envelope masses < 0.2 Mq (Fig 1 1 0b and having a compact star 
companion. The compact nature of their companions allows for ex¬ 
treme stripping in a tight binaries which is not possible otherwise. 
We find several cases of Case BB RLO with P or b,i < 2 days lead¬ 
ing to ultra-stripped SNe, i.e. exploding stars which are almost pure 
naked metal cores with only a tiny amount of helium envelope left, 
in some cases only ~ 0.005 Mq for EC SNe and ~ 0.02 Mq 
for Fe CCSNe. In general, we find that the progenitors of EC SNe 
are stripped slightly further than those of Fe CCSNe, meaning that 
our EC SNe are expected to be mainly Type Ic SNe, whereas the 
Fe CCSNe from these systems can be either of Type Ic or Type lb. 
The amount of helium remaining in the envelope is found to be 
correlated with the initial orbital period (Figs. l9landl 16b. and thus 
also the final pre-SN orbital period, of the binary system. In addi¬ 
tion, there is a correlation of remaining total envelope mass with 
the final core mass of the exploding star (Fig. m. 

From our modelling we expect ultra-stripped SNe to pro¬ 
duce NSs (the second formed NS in DNS systems) within a rel¬ 
atively wide range of masses (possibly even in the entire range 
of 1.10 — 1.80 Mq, Fig. QT}. Besides forming low-mass NSs of 
1.10 — 1.40 Mq from almost naked metal cores barely above the 
Chandrasekhar mass, we find cases where we expect ultra-stripped 
SNe to produce more massive NSs in tight orbits, originating from 
metal cores of evolved helium stars with initial masses of 4—6 Mq. 

On the basis of our modelling, we conclude that EC SNe 
only occur for a narrow interval of initial helium star masses of 
Mne,i = 2.60 — 2.95 Mq depending on P orb ,i (Fig. [2j. The gen¬ 
eral outcome of our investigated binaries, however, is an Fe CCSN 
above this mass interval and an ONeMg WD or a CO WD for 
smaller masses. We discuss one peculiar model where a 2.60 Mq 
helium star donor completely shuts off core nuclear burning as a 
consequence of extreme mass transfer (Figs. \m and we followed 
its subsequent evolution as a CO WD to the UCXB-phase where it 
filled its Roche lobe with an orbital period of about 30 seconds. 

The maximum mass-transfer rate during stable RLO from all 
the binary helium stars is close to ~ 10” 4 Mq yr -1 within a factor 
of a few (cf. Figs. EHH. The duration of the mass-transfer phase 
is anti-correlated with the initial P orb , which results in the predic¬ 
tion of a correlation between P or b and P sp i n for the mildly recycled 
pulsars in DNS systems. The reason is that in closer orbits the du¬ 
ration of the mass-transfer phase is longer (Figs. [T3landl 16b. which 
results in more mass being accreted by the NS, leading to more ef¬ 
ficient recycling and thus a smaller value of P sp i n (and assuming 
this relation to survive the dynamical effects of the SN explosion). 

The main channel for producing DNS systems which will 
merge within a Hubble time is via ultra-stripped SNe. It will 
therefore be useful to compare upcoming DNS merger rates from 
LIGO/VIRGO with future observational limits on the detection rate 
of ultra-stripped SNe. The result of this comparison may help to 
answer whether, and to which extent, ultra-stripped SNe also oc¬ 
cur in binaries with a WD or a black hole companion star. How¬ 
ever, this step requires a clear method to distinguish the light 
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curves of ultra-stripped SNe from other fast decaying and relatively 
weak li ght curves, e .g. resulting from partial explosions (' .la’ SNe, 
iBildsten et al.ll20 07h or the accretion-induced collapse (AIC) of a 
WD l iDessart et al .2006b - see also Section l4~2l for discussions and 
references to spectroscopic signatures of diverse SNe. 

Ultra-stripped SNe have small amounts of ejecta mass, often 
with a low binding energy (Fig. d. We therefore speculate that 
ultra-stripped SNe (both EC SNe and, at least, Fe CCSNe with rel¬ 
atively small iron cores) lead to fast explosions and often produce 
small kicks. 

We also performed a brief investigation of helium star-NS bi¬ 
naries with P 0 rb,i > 10 days (Fig. d leading to the formation 
of wide-orbit DNS systems after the second SN explosion. These 
DNS systems often have P OT b > 20 days and are characterized by 
inefficient recycling, often leading to spin periods of the marginally 
recycled pulsar in excess of 100 ms. Furthermore, their eccentrici¬ 
ties are expected to be large (0.2-0.5, even for symmetric SNe with 
no kicks) because such exploding stars eject fairly thick envelopes 
and their pre-SN orbits are relatively wide, and thus weakly bound. 
In a future publication (Tauris et al., in prep.) we analyse in more 
detail the observed distributions of orbital periods, spin periods, 
eccentricities and systemic velocities of DNS systems and confront 
them with simulations and expectations from ultra-stripped SNe. 

In a few cases, envelope expansion triggered by a helium shell 
flash following ignition of oxygen burning leads to a CE-like sit¬ 
uation with the NS embedded in the outer helium envelope of the 
donor star. However, we demonstrate that such events will not lead 
to a runaway mass-transfer (and a dynamical instability) which 
would cause the NS to merge with the core of its companion star. 
The reason is that the gas in the expanded envelope is so dilute that 
the timescale of spiral-in (cf. Appendix) often exceeds the remain¬ 
ing lifetime of the pre-SN star by a factor of (at least) 100. In addi¬ 
tion, the binding energy of this expanded envelope is much smaller 
than the release of orbital energy from a slight spiral-in in these 
tight orbits. Thus we conclude that it is safe to ignore the expand¬ 
ing envelope at the very late stages of evolution when determining 
the final outcome of such models. 
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APPENDIX A: TIMESCALE OF NEUTRON STAR 
SPIRAL-IN CAUSED BY A HELIUM FLASH 


For the purpose of determining the fate of binaries in which the 
donor star experiences a vigorous helium shell flash after (or close 
to) oxygen burning has been ignited, we briefly outline an analyt¬ 
ical investigation to demonstrate that the timescale of spiral-in is 
much longer than remaining lifetime of the pre-SN star. 

As an example, we shall analyze the situation for the system 
with initial parameters: Mue,i = 2.9 Mg, P OI b,i = 0.1 days and 
Mns = 1.35 Mg . This systems was studied in detail in Paper I 
dTauris et all2013D . The vigorous shell flash gave rise to numerical 
problems for our stellar code and we alleviated this problem by 
evolving the donor star as an isolated star during the last few years 
of evolutior@. Using this method, we could follow the expansion 
of the donor star during the flash. The donor star expanded from 
R = 0.40 to 2.41 Rq, which would engulf the NS inside the dilute 
helium-rich envelope of the donor. As we briefly argued in Paper I, 
this star would undergo an Fe CCSN within only 10 yr and thus we 
could neglect the spiral-in caused by the frictional drag acting on 
the NS prior to the explosion. Here we elabor ate a bi t more on this 
argument and follow the method suggested bv lDewi & Pod ( 2003 ) 
to evaluate the timescale of the spiral-in of the NS. 

The loss rate of orbital energy, |J7 or b| as a result of the drag 
force, Tdrag acting on the NS inside the dilute envelo pe of the 
evolved helium s tar is roughly given by iBondi & Hovlel fi 9441 : 
Shima et al ]fl98l) : 

|J5 or b| — -^drag V — A -Pdrag V — Up) 7tRa.cc P V , (Al) 


where the NS accretion radius is: 


R 


acc 


2GMns 

V 2 + Vs 


(A2) 


Here G is the gravitational constant, v is the relative velocity be¬ 
tween the NS and the CE, v 3 = y/'y P/p is the speed of sound, 7 
being the adiabatic index (5/3 for a monatomic gas), P and p are 
the local pressure and mass density provided by our stellar mod¬ 
els, cf. Fig. AL Ut 1 ) * s a scaling factor determining the dissipa¬ 
tion rate of E OI v due to the drag. Here we chose Up) = 2.5 
for the supersonic case (in our models v is larger than v s by a 
factor of 3 — 10). Similarly, the rate at which the orbital energy, 
E OI b = — GMH e MNs/(2a) is dissipated is given by: 


| -E/orb | 


G Afee Mns . 

---a 


(A3) 


where a is the orbital separation between the NS and the centre 
of the evolved helium star. Assuming Keplerian motion for the in- 
spiralling NS (and neglecting for a moment corotation of the CE): 
v(a) = yjG (MHe + Mns)/o. and this enables us to equate the 
expressions in equations ED and El- 

Integration then yields the timescale of the spiral-in, for ex¬ 
ample, from the current location of the NS, a 0 until it reaches the 
metal core boundary of the evolved helium star at radius, R core : 


id 


ecay 



_ G Mue Mns _ 

2 a 2 £(/x) 7 r R acc (a ) 2 p(a) v(a) 3 


da. 


(A4) 


2 Note, there is a typo in the second paragraph of Section 2 in Paper I. 
The off-centre oxygen burning was (re)ignited at t = 1.854 353 Myr (not 
t = 1.854 553 Myr), i.e. only some 16 yr in evolutionary time after we 
took the star out of the binary and evolved it as an isolated star to investigate 
its expansion due to the flash. When t had increased by an additional few 
years (3-4 yr) the computations finally broke down at t = 1.854 356 Myr. 
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Radius coordinate, r (R Q ) 

Figure Al. The mass-density profile of the evolved 2.9 M& helium star 
(Porb. j = 0.1 days and = 1-35 Mq) investigated in lTauris et alj 
i2013l) , see also Table [7] for further information on this system. The total 
mass of the expanding outer envelope is 2.58 X 10 -3 Mq and its average 
mass density is 2.6 X 10 — 4 gem -3 . The location of the NS is shown. 


Radios coordinate, r (R„) 

Figure A2. The orbital decay time, tdecay as a function radius coordinate. 
Blue curve: tdecay is calculated assuming that the NS spirals in all the way 
to the centre (r = 0) of the evolved helium star. Red curve: NS in-spiral 
until the boundary of the metal core (r = 0.024 Hr./). Orange curve: the 
local characteristic orbital decay timescale, r = E OI p / 7r or i, . See text. 


In Fig. A2, we have plotted tdecay for the system discussed in Pa¬ 
per I. The blue line corresponds to R co re = 0 (i.e. the spiral-in 
timescale for complete merger to the centre of the star). The red line 
(similar to the blue line for the majority of the computation) shows 
the timescale for spiral-in until reaching the boundary of the metal 
core, .Rcore = 0.024 R@. The orange line simply shows a local 
characteristic orbital decay timescale defined as: r = E or b/E or b. 
One can see that for the present location of the NS, r « 3500 yr 
(3.8 x 10 48 erg/3.2 x 10 37 ergs -1 ). However, r is underestimat¬ 
ing the rate of spiral-in, because the density increases inward, and 
therefore r is an overestimate of the more realistic orbital decay 
time. This is seen in comparison with the blue line which yields 
fdecay = 600 yr for the location of the NS in our final model 
(ctf = 1.02 Rq and P or b,{ = 0.070 days). 

During the helium shell flash it is only the outermost part of 
the star which expands. In the specific model discussed here it is the 
outer 2.58 x 10 -3 Mq of material. This outer envelope is expanded 
by a factor 10 in radius during the flash (reaching 2.41 Rq) and 
explains the relatively weak drag on the NS due to the dilute gas. 

To conclude, it is safe to assume that t decay ~S> t collapse, where 
fcollapse is the remaining lifetime of the evolved helium star un¬ 
til the SN explosion. Therefore, our calculated models will lead to 
ultra-stripped stars exploding without the NS companion being part 
of their core regiorf^. This conclusion is further strengthened if one 
includes the deposition of orbital angular momentum into the enve¬ 
lope which wou l d cause it to cor otate and thus reduce the drag force 
on the NS jLivio & Sokerfl988l) . leading to larger values of tdecay- 
In addition, one must bear in mind that in the example demonstrated 
here, we considered a very conservative case with fcoiiapse — 10 yr 
(which is still small compared to tdecay = 600 yr). However, in 
other models we even have t C oiia P ae ~ a few days, cf. Section [3Tl 

Even in case the NS would experience a slight drag from 
spiral-in, the released orbital binding energy in this tight system 
is more than sufficient to effectively remove the loosely bound en¬ 
velope (Ai? or b Ebind, by more than two orders of magnitude). 


3 In case the NS would have merged with the core of the pre-SN star it 
would presumably undergo hypercritical accr etion , leading to formation of 
a black hole and possibly a 7 -burst like event toermer & Atovanll2006i) . 
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